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SU-IMARY 


The objective o f this study wc»s to formulate a conceptual design for 
an ion thruster that would meet the anticipated performance requirements 
of large space systems. The principal thruster performance goals for the 
design analysis performed under this contract were a thrust of 0.5 N, 
operation on inert gas propellant (xenon or argon), operation at a 
specific impulse of 3530 sec for xenon and 6076 for argon, and achieve- 
ment of a total thruster efficiency of 74% for xenon and 68% for argon. 
The design study was totally analytic, drawing heavily on the analytic 
models for inert gas thruster operation that have been developed at 
Colorado State University. As a consequence, the thruster discharge 
chamber makes use of a multipole magnetic confinement configuration for 
which scaling and analytic description have been demonstrated. To 
achieve the thruster efficiency goals, a beam diameter of approximately 
50 cm was found to be necessary. Having determined that this was the 
major design driver, tractable solutions to all of the remaining design 
problems were worked out (in concept). The resultant conceptual design 
for the "0-cm inert-gas thruster provides a promising basis for further 
engineering and confirmation of untested design approaches (such as 
multiple hollow cathode operation). 
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SECTION 1 


INTRODUCTION 


Ion thruster development has concentrated primarily on mercury 
electron bombardment thrusters that accelerate mercury ions electro- 
statically. The technology for thrusters producing 8-cm and 30-cm- 
diameter ion beams (producing thrust levels of 5 mN and 125 mN, respec- 
tively) has reached a relatively advanced status. These thrusters are 
now considered ready for application as auxiliary propulsion (8-cm 
thrusters) for satellites in earth-orbit and for primary propulsion 
(30-cm thrusters) for interplanetary exploration. Projecting to the 
future, large space systems (LSS) are contemplated for near-earth space 
platforms (earth orbital systems) having a variety of functions. To be 
cost effective, ion thrusters would have to produce greater thrust (per 
thruster), have long life, and be compatible with other LSS functions. 

The work described in this report represents an initial analysis and 
conceptual design for an ion thruster to meet the projected requirements 
of LSS. 

Environmental concerns over the relatively large quantities of 
propellant that would be required for an LSS preclude using mercury as a 
propellant, and argon and xenon were specified as the required alterna- 
tive propellants. This represented the major design constraint, and the 
design analyses were performed with the objective of showing the capa- 
bility of the thruster concept developed for meeting the performance and 
operational characteristics specified as the contract goals. At the 
outset, these goals (described in Section 2) seemed to be a challenging, 
if not overly ambitious, extrapolation from the state of the art. At 
the completion of the analyses, however, the contractual goals appear 
realizable with some reasonable and realistic assumptions about scaling 
the thruster design using experimentally verified principles and cor- 
relation parameters. With regard to performance analyses and experimental 
correlations, the work reported herein draws heavily on research conducted 
at Colorado State University^ ^ (unde.. NASA grants). 
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SECTION 2 


DESIGN ANALYSES 


The analysis of ion chruster performance must take into a -count 
several figures of merit; (such as efficiency, operating range, and 
reliability) that are interdependent. Consequently, we began by evalu- 
ating the design goals with respect to which characteristic was most 
difficult to satisfy. This section describes the design goals, identi- 
fies the design drivers, and outlines the computational procedures used 
to obtain the required thruster characteristics. 

A. DES'.'GN GOALS 

The design goals specified in the statement of work for this contract 
are summarized in Table 1. The dominant independent parameters specified 
in this table are thrust, specific impulse, and type of propexlant. As 
shown below, all the other quantities have to be derived from these three 
quantities using analytic models for thruster performance. Since analytic 
models for ion thruster performance are only reliable for certain thruster 
configurations, the design features specifying module scalability, 
variable thrust, and specific impulse are more restrictive than might 
appear. Only the multipole discharge chamber geometry described by 

fi 

Isaacson has been scaled analytically and verified experimentally. Scala- 
bility is the major design driver for discharge chamber configuration. 

Other design drivers are discussed in the following section. 


B. DESIGN DRIVERS 

We initially thought that thruster lifetime would be the major 
design driver dictating thruster diameter. But computation of thruster 
efficiency (based on specifying the design parameters to satisfy screen 
grid lifetime) showed that thruster efficiency goals could not be satis- 
fied using this approach; therefore, satisfying efficiency became the 
major design driver. Efficiency computations by Kaufman (described in 
Appendix A and also in Ref. 12) were useful in defining the thruster 


13 


Table 1. Specifications for Design of a Thruster for 
LSS Application (Design Goals) 


Propellant: inert gases argon or xenon 

Thrust: 0.5 N 

Specific impulse: 

Argon, 6076 sec 
Xenon, 3530 sec 
Thruster efficiency: 

Argon, 68% 

Xenon, 74% 

Lifetime: 15,000 hr 

Design features 

• Specific impulse and thrust variable ±50% 

• Module scalable from half to double nominal 

• Constructed from commercially available materials 

• Concepts ultimately applicable to space flight 

• Unitized construction for oossible replacement of critical 
components 
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diameter through the relationships shown in Figures 1 and 2. These 
computations were made under the assumption that the thruster can be 
operated at a discharge voltage sufficiently low to neglect the effects 
of doubly or triply charged ions. Unless the ratio of the length of the 
discharge chamber to its diameter is very small, multiply charged ion 
effects are found to be appreciable for operation with xenon as the 
propellant. Accounting for these effects was facilitated using the 
analyses and design curves prepared by Wilbur and presented in Appendix B. 
Using overall thruster efficiency as the major design driver, the com- 
putational technique used in defining the design is described in 
Section 2.C. 

C. COMPUTATIONAL PROCEDURES 

To compute the theoretical performance of an ion thruster, the 
quantities identified in Table 2 must be specified or computed using the 
relationships derived by <aufman and Wilbur^ and described in detail in 
Appendices A and B. Having noted from the cross-plots of Kaufman's 
results shown in Figures 1 and 2 that the beam diameter for maximum 
thruster efficiency is ^44 cm for the values of thrust and specific 
impulse under consideration here, it remained to determine the optimum 
discharge chamber length. For computational purposes, we assumed that 
the cyclic pole-piece/anode dimension is 1.1 cm (the same as that used 
by Isaacscn at Colorado State University). The thruster operating 
parameters and total efficiency were computed as a function of discharge 
chamber length to obtain the values listed in Table 3; the procedure is 
illustrated by the sample calculation for xenon propellant which follows. 

A nominal beam diameter of 50 cm has been assumed for these computations. 
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Figure 1. Thrust versus bean diameter for argon propel- 
lant with specific impulse as a parameter 
(thruster operated for maximum efficiency). 
Dashed line is specific impulse design goal. 
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Figure 2. Thrust versus beam diameter for xenon pro- 
pellant with specific impulse as a parameter 
(thruster operated for maximum efficiency). 
Dashed line is specific impulse design goal. 
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Table 2. Definition of Syebols Used for the Quantities 
Necessary for Computing Thruster Performance 
Characteristics (Including Assigned Values 
where Applicable) 


Quantity 

Symbol 

Volu* 

Propellent aasa 

a 

o 


Specific lapulu 

'.p 


Thrust 

t 


IMS disaster 

s 


team area 

s 


Discharge cheater length 



Accelerator grid transmission 


0.2 

Screen grid transmission 

♦. 

o.r 

Screen grid thickness 

t 

s 

0.0S1 cm 

lean voltage 

v b 


team current 

J b 


team power 

% 


Bean unlforatty paraaste- 

r b 

0.»5 

Total thrust correction (actor 

r 


Discharge voltage 

v d 


Emission current 

J E 


Discharge power 

p o 


Neutral propellant flow 

^O 


Propellant efficiency (Manured) 



Electrical efficiency 

n . 


Thrust correct ion factor for Heaa divergence 

r l 

0.0M 

Thrust efficiency 

'V 


Thruster (screen grid) llfetiM 

z 


Primer y electron region surface area 

S 


Prlaer slsctron region voluM 

V p 


E con teaperature (Maxwellian) 

T « 

4 oV 

Plasma uniformity factor 


5 or 7 

Neutral atoa velocity 

y o 

•aoo’K 

Doha velocity (Iona) 

»♦ 


Neutral loss pa ream ter 

P m. 

1.4 Ur); 0.41 (So) 

Discharge loss paraaetar 

p ot 

47 Ur); 42.2 do) 

Ionisation efficiency 

n i 


Pius of singly charged Ions 

r * 


Flux of doubly charged ions 

r ~ 


Discharge plasaa paraaetar 

e ) 


Ionisation rata factors 

*i 


Ratio of doubly to singly charged ion flux 



Thrust correction factor for doubly charged ions 

■- 


Primary electron confinement paraaetar 

m* 

54 n 10** Toola-m 

Anode area 

A A 


Maximum anode current 

j am 


Specific discharge power 

P St> 


Minimum discharge voltage 

V 0m 


Sputtering yield for singly charged Ions 



Sputtering yield for doubly charged Ions 



Density of screen grid Mterlal 

p 

4.4 a 10 22 cm" 2 

Total acceleration voltage 

mm 


tatio of beaa voltage to total acceleration voltage 

EM 








Theoretical Operating and Performance Quantities Computed as a Function of 
Discharge Chamber Length for a 50-cm (beam diameter) Thruster 
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1. Compute areas and volumes using: 


Bean diameter 


d b 


50 cm 


Chamber length = i . * 8.1 cm 


a. Primary electron region surface area 

2 \ /% 

I 4- TtA 9 * iM I 4- L 

d 


A = 2 
P 


(■v) + "Vd * ” d b(r + 0 ' °- 5182 " 2 


b . Anode area 

"djj 2 

A A = — + 11 Vd - A cath = 0-3154 


( 1 ) 


( 2 ) 


area of cathode plus other vacant area 


Beam area 

.2 


TTd 




b 2 

- = 0.1963 m 


(3) 


d. Primary electron region volume 

wd?fc. T 

Y = — — = 0.0157 m 3 
P 4 


(4) 


e. Volume-to-surface area ratio — primary electron region 


!e._ 

A 2d. 
P 


Vd 


b + 4t d 


= 0.0303 m 


(5) 
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n* ■ 


o r:c: 

W fi. 


i 


Compute neutral propellant flow rate 


I 


Te 

m g 
sp o o 


10.6 A 


» 


-19 

where T « thrust (in N), e * 1.6 x 10 C, I = specific 

impulse (in sec), m - propellant atom mass (In kg), and 

g ■ 9.8 m/sec^. 
o 


Compute ionization efficiency (defined here as the ratio 
of the ion loss rate to the neutral atom supply) 


p niM 


n ‘ ‘ 1 ' vyy 


= 0.95 


where P NL is the neutral loss parameter (from Ref. 6), and 
is the transparency of the accelerator. 


Compute the plasma parameter "C" for F + = 5 



Compute chamber specific discharge loss level 

P SD ‘ P DL £ - 165 eV/i °" ' 
b 

where P = 62.2 eV/ion (obtained from Ref. 6 and adjusted 
for 4> s = 0.7). 
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(7) 


( 8 ) 
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6. The minimum discharge voltage is determined by the 

criterion for multipole thruster anode area, which determines 

the maximum anode current J A11 by the relation 

AM 


kT n A a 
e e A 

AM " 13/Bdx 


( 10 ) 


coupled with the definition of specific discharge power 


= ? SD J b 

0111 " Wb 


(ID 


and the expression for the discharge chamber electron 

•k 

density (n^) 


j b pr 

n e - n i - eA b 4, s V kT e 

where is the electron temperature, k is Boltzmann's 

constant, A^ is the anode area, J A is the anode current, 

and J, is the beam current. Combining the last three 
b 

equations yields the expression for minimum discharge 
voltage determined by the anode current limitation: 


( 12 ) 



24.8 


•k 

The correct expression for density in this case is 


(13) 


n 

e 


J b . 

f—^o 

i+^W 

’ e Vs V 

kT e ( 1 + ”p /n n ) 

i + 2 (r ++ /r + ) 


but the equation given in this case should be an adequate approximation. 
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7. Having determined the minimum value for discharge voltage, 

one can estimate V and V -V such that 
U U K 


V D > 24.8 


( 14 ) 


For the proposed design, since V -V„ was chosen to be 28 V, 
and was assumed to be 6 V, it follows that V = 34 V. The 
remaining design values are found from the desigu curves 
shown in Appendix B. Using the value of 0.53 computed for 
the plasma parameter, it is found that 

1^71^=0.27 ; (15) 


from Figures B-l and B-3, it is found that 
V, r\[ = 686 

D 1 

= 760 . 

b 

Figure B-4 produces 

— = 12.9 
n i 

J. '« 12.3 , 

b 


(16) 

(17) 


(18) 

(19) 


assuming no beam divergence. Including the effect of 
beam divergence (F t = 0.985), 


J L =12.5 . (20) 

b 

From Figure B-5, the thrust correction factor for doubly 
charged ions, a, is found to be 

a = 0.892 . (21) 
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The total thrust correction factor becomes 


Y ■ a F ■ 0. 879 , 


the beam power is 


p b * v b J b ' 9500 • 


the discharge power is 


P D * Ve * J b P SD * 2062 • 


and the total power becomes 


P TOT P b + P D 11 » 563 


The electrical efficiency is 


% ' V P D 


= 0.821 


The propellant efficiency that would be measured in an 


thruster is 


n = ~ - 1-18 . 

m J- 
o 


The overall thruster efficiency is 


n T = Y = 0. 747 

T e tn 


( 22 ) 

(23) 

(24) 

(25) 

(26) 

operating 

(27) 

(28) 


This procedure was repeated for all of the values of listed in 
Table 3. 


OUi'cUiWu t to 

OF POOR QUALITY 


If the losses associated with a neutralizer are included, 

then J has to be increased by the neutralizer flow 
o 


J 0T " J o + J 0N 


= 10.6 + 1.0 = 11.6 


* J b 
0 _ = 


m J 


= 1.08 


OT 


Electrical efficiency is modified to be 


p +p + p = 0*806 

D b N 


P. T = 16(J, + 1.5) = 224 
N b 


The overall corrected thruster efficiency becomes 


( 29 ) 


(30) 


(31) 

(32) 


n T = 0* 0* Y 2 = 0.672 . (33) 

T me 

Note that the values of n shown in Table 3 for operation with Xe 

m 

propellant are always greater than 100% because of the relatively high 
percentage of doubly charged ions. Surprisingly, as shown in Section 2.D, 
this high percentage of doubly charged ions does not present a screen- 
grid lifetime problem. On the other hand, the relatively high discharge 
voltage required for operation with argon (to satisfy the maximum anode 
current condition) does present a severe lifetime problem. Triply 
charged ions could be appt 3ciable if the discharge voltage were greater 
than 34 V for Xe operation, and if the ionization efficiency (n ) for Ar 
operation were to exceed 0.9. Neither of these conditions occurs at 
the design operating point selected here and therefore the effect of 
triply charged ions can be neglected. 
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A plot of overall thruster efficiency (n T ) as a function of discharge 
chamber length (8,^) , given in Figure 3, shows that the best performance 
(efficiency) is obtained at a different value of 8.^ for operation with Ar 
than for operation with Xe. Specification of an of 16.2 cm produces 
performance that meets the contractual goals for operation with both 
Xe and Ar. 



0 5 10 15 20 

DISCHARGE CHAMBER LENGTH, cm 


Figure 3. Comparison of thruster efficiency (corrected 
for doubly charged ions and beam divergence) 
as a function of discharge chamber length 
for operation with xenon and argon propellants. 
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D. THRUSTER LIFETIME COMPUTATION 


0°"' 

c; i 


n f ;■ r -i 

r-H -v 


'J\-U 


The principal life-limiting factor has not yet been determined for 
inert gas ion thrusters; however, it is reasonable to assume that compo- 
nent development will result in the same situation for inert-gas thrus- 
ters as for the mercury ion thruster. The screen grid electrode of 
the ion optics assembly must be fabricated from relatively thin material 
to achieve efficient ion extraction. If the screen grid is sufficiently 
thin, it is thought that ions may, in fact, miss the electrode, and hence 
ion sputtering would not present a lifetime problem. Until this suppo- 
sition has been verified, the thinning of the ion optics screen grid 
electrode by ion sputtering is considered to be the principal wearout 
mechanism. The expression used to calculate the useful lifetime, t, is 


T 


,°.9 t a e Ve F b pW i ♦ r++ /r + \ 

\ 3t00 Vo As + + s+A nj' 


(34) 


where T is the grid lifetime (in hr), t^ is the grid thickness (in cm), 

2 S 

is the area of the beam (in cm ) , ^ is the screen grid open area 

fraction, is the ion beam flatness parameter, p is the density of 

the grid material, S + is the sputter yield of the grid material corres- 

| 1 1 ( 

ponding to the energy of singly charged ions (%V^), and S is the 
sputter yield for doubly charged ions (^2 V^). The factor 0.9 in the 
numerator corresponds to an assumed condition of screen grid failure 
when 90% of the initial grid thickness has been eroded away. Using a 
grid thickness of 0.051 cm, a screen grid open area fraction of 70%, a 
beam flatness parameter of 0.95, and other parameters defined previously, 
F.q . (34) becomes 


'Jh 

A, 


8.52 x 10 


-3 


I 1 + r +f/r+ \ 
\s + + s 44 r++ / r+ / * 


(35) 
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As shown in Appendix B, F /F, is a function of the plasma parameter C 

+ 

and the electron energy (V - V„). The sputtering coefficients S and 
++ UK 

S are functions of V, (see Table 4). If V_ =* (V_ - V„), then the life- 

b L) UK. 

time parameter can be computed as a function of the plasma param- 

eter, as shown in Figures 4 and 5. If V D >> (V^ - V^) , then it is nec- 
essary to use the pertinent values of S and F | , /F ( . respectively, in 
Gq. 35. 

If a thruster is operated on xenon propellant at values of 
(Vp - V^) > 32 V and relatively higher values of n^, then the grid life- 
time could be altered appreciably by triply charged ions. The expression 
for the grid lifetime is 

8.52 x 10* 3 A 

D 

I * 


i + r -H-/r+ + r +++ /r + 


(36) 


s + s r. ./r. + s r , , , /r , 


Table 4. Sputtering Coefficients for Sputter Yields 
of Singly and Doubly Charged Ions on Molybdenum 


as Functions of Discharge Voltage, 


Gas 

V D 

S + x 10 -6 
atoras/ion 

S ++ x 10" 3 
atoms/ion 

Ar 

40 

40 

40 


35 

20 

16 


34 

15 

15 


32 

9 

9 

Xe 

35 

9 

5. 3 


34 

7 

4 


32 

5 

3 


30 

2 

1.5 


28 

1 

0.5 
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GRID LIFETIME PARAMETER 







GRIO LIFETIME PARAMETER 






cf poor. 


q’.j hv : c ( 


+ 44~ 4 -f-f 

The quantities S , S , rnd S are the sputtering yields for singly, 

doubly, and triply charged ions, and A, is the beam cross-sectional area. 

b 2 

For a beam cross-sectional area of 2000 cm and an ionization efficiency 
r) of 0.96, the grid lifetimes are given in the following table along 
with sputter yields for the triply charged and the corresponding lifetimes 
fjr the case where triply charged ions are not considered: 


E, eV 

,4-H- 

a 

b , 

S 

t , hr 

t , hr 

35 

0.09 

11,400 

12,500 

32 

0.045 

25,000 

25,000 


With triply charged ions 
^Without triply charged ions 


The above table shows that triply charged xenon ions can have a small 
effect on grid lifetime when the discharge-keeper voltage difference is 
35 V. This effect is eliminated when this voltage is dropped to 32 
or less. 
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E. 


MAGNET J FIELD CALCULATIONS 


It is anticipated that the thruster body, including the magnets and 
pole-pieces, will have to operate at a relatively high equilibrium 
temperature. Consequently, the magnet material selected should have 
superior resistance to demagnetization. Of the magnetic materials 
available, ALN1C0 8 appears to have the best characteristics with regard 
to demagnetization. In comparison with ALN1C0 5, the material used for 
30-cm mercury ion thrusters, ALN1C0 8 has a better resistance to demag- 
netization. For example, in reducing the magnet length-to-diameter ratio 
(L/D) from infinity to a value ,»f 4, ALNICO 8 loses only 15% flux, whereas 
ALNICO 5 loses 25%. Figure 6 shows a demagnetization curve that is 
representative of ALNICO 8 material. Assuming a value of 4 for the 
L/D ratio, a flux density of 0.75 T is considered reasonable. This value 
will not be sensitive to small changes in L/D. 

In a discharge chamber that employs a multipole magnetic confinement 
configuration, the magnetic circuit is composed of ring-shaped pole-pieces 
for the wide-wall discharge chamber boundary with short permanent magnets 
spaced between the »-tngs, as shown schematically in Figure 7. The upstream 
end of the discharge chamber is composed of concentric cylinder pole- 
pieces with the permanent magnets spaced around each annulus. The pole- 
piece thickness is •‘equired to be approximately 1/4 of the magnet diameter 
to avoid saturation where the magnet contacts the pole-piece. The milled 
recess that holds the magnet in position does not change this requirement 
because it Is the thickness of the material that surrounds the end of the 
magnet that is important. It is important that magnets from adjacent 
"cusps" are not aligned on opposite sides of the pole-pieces. Although 
this is relatively easy to satisfy on the side walls, magnet placement 
in the upstream end pole-pieces requires some compromise between achieving 
a uniform spacing between magnets and separating the point of contact of 
magnets on opposite sides of the pole-pieces. 

The objective of the computations performed here is to determine 
toe number of magnets required to generate the desired confinement fields. 
The procedure for calculating the magnet requirements has already been 
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developed at Colorado State University and reported in Ref. 4 and 11. 

As discussed in Ref. 11, the curvature of the side wall pole-pieces is 
most important. But since the curvature effects for a 30-cm-diameter 
chamber were found to be small, these effects should be even smaller for 
the proposed 50-cm chamber. From Ref. 4 (30-cm chamber), the curve 
describing the magnetic field variation in the pole-piece region can be 
integrated to show that about 20% of the flux falls within the inner 
surface of the anode (see Figure 7), with this inner surface assumed to 
have the same diameter as the inner edge of the pole-piece (as shown 
in Figure 7). For the proposed 50-cm design, the radial depth of the 
pole-piece is 3.5 cm, while the axial spacing is 2.7 cm. This depth-to- 
spacing ratio would be expected to reduce the flux inside the anode 
inner edge to about 16% of the total, compared to the 20% for the 2.5-cm- 
deep design of Ref. 4. 

For adequate primary electron containment, the flux integral over 
an anode (toward the axis) should be about 1.5 times the theoretical 
value given by 


Bdx = 6.74 x 10~ 6 E 1/2 , 

e 

where E^, the primary electron energy, is in eV. For a maximum primary 

electron energy of 35 eV, the theoretical integral should be about 

40 x 10 ^ T-m (40 G-cm). Something around 60 x 10 ^ T-m should therefore 

be used in the design. For the low discharge voltages of interest, a 

low value of this integral is desired to maximize the electron diffusion 

limit. An integral of 54 x 10 ^ T-m (1.4 times the theoretical value) 

has therefore been selected. The total flux at about a 50-cm diameter 

(above the anodes) is thus tt ; 0.5 x 54 x 10 1 Wb (neglects curvature 

effect). With this flux, 16% of the total, the total flux should be 

about 5.3 x 10 Wb. The flux through one 6.35 ram (1/4 in.) diameter 

- 3 2 

magnet at the assumed flux density of 0.75 T is (6.35 x 10 ) (tt / 4)0. 75, 
or 2.375 x 10 Wb. The number of magnets for one pair of pole-pieces 
is therefore 5.3 x 10 ^ / 2 . 375 x 10 ^ = 22.3. One should not be 


35 


overimpressed with the accuracy of this calculation. A value of 22, 

though, will be used here for further calculations. 

For the upstream end of the discharge chamber, the number of 

4 

magnets must be varied with the anode diameter. Note that the inner 
diameter of the side wall anodes is used, not the mean diameter. The 
justification for this choice is simply agreement with the experimental 
fields observed. Using this means of varying magnet number, one finds: 


Anode Number of 

Diameter Magnets 


Side-wall anodes 

49.7 cm 

22 


End-wall 

anode 

47.0 

20.8 

+ 21 

2nd -wall 

anode 

41 . 6 

18.4 

+ 18 

End -wall 

anode 

36.2 

16.0 

+ 16 

End-wall 

anode 

30.8 

13.6 

+ 14 

End-wall 

anode 

25.4 

11.2 

+ 11 

End -wall 

anode 

20.0 

8.9 

+ 9 

End -wall 

anode 

14.6 

6.46 

+ 6-7 

End-wall 

anode 

9.2 

4.1 

+ 4 


It may be necessary to change the diameters of the upstream pole-pieces 
slightly, but the method used to select magnet numbers should be clear. 
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F. 


THERMAL AND MECHANICAL DESIGN CONSIDERATIONS 


The primary concern in the mechanical design for operation at the 
beam current levels being considered is the rejection of the thermal 
energy from the discharge and the relatively high operating temperatures 
that will have to be accommodated by the thruster structure and ion 
optics assemblies. Our initial analyses have been based on the concept 
that heat rejection can be obtained by direct radiation from the compo- 
nents heated and that heat pipes or other heat transfer devices will not 
be required. 

In the LeRC/Hughes J-series 30-cm thruster, the thruster structural 
elements act as a heat shield between the anode and the ground screen, 
which would require extremely high anode temperatures. A multipole 
magnetic field geometry is favored for meeting the design goals of this 
study, and the possibility exists that the discharge chamber can be 
constructed with a configuration that permits the multipole anodes to 
radiate directly through the ground screen. Figure 8 shows several 
typical pole-piece/anode sections. The anodes are formed to have a "T" 
cross section; the leg projects into the ionization region, and r he 
"cross" of the "T" forms the discharge chamber wall. Electrical isola- 
tion of the interior pole-piece (Figure 8) from the discharge chamber 
(which is at anode potential) is obtained by providing a gap between the 
interior pole-piece and the exterior pole-niece that holds the magnets. 

Previous thermal analysis of LeRC/Hughes 30-cm mercury ion thrusters 
was based on the assumption that approximately one-half of the discharge 
power is deposited in (and has to be radiated from) the discharge anode. 
In a multipole thruster, anode currents should be relatively uniformly 
distributed around the internal surfaces of the discharge chamber 
(predominantly to anodes) and consequently it is reasonable to assume 
that the discharge power will be uniformly radiated from the anode 
surfaces. The worst-case condition occurs for operation with argon 
because of the higher beam current and discharge power required. Pre- 
vious thermal modeling ot the 30-cm mercury ion thruster apportioned 


37 



ORIGINAL PAGE IS 
OF POOR QUALITY 


9573-40 



Figure 8. Typical magnet/pole-piece cross section showing anode 
configuration that provides for direct radiation of 
thermal input. 
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the power inputs to the discharge chamber components (approximately) 
as follows: 

Anode 105 

Screen 20 J, 

b 

Backplate 20 

Lost through 40 

optics 

185 J, 
b 

If the specific discharge power is apportioned accordingly for the 50-cm 
multipole thruster, the power input to the components will be: 

Anodes 140 J, 

b 

Screen 33 J, 

b 

Lost through 33 

optics 

206 J u 
b 

In the multipole thruster, the entire interior surface of the discharge 

chamber is effectively anode area. For the 50-cm multipole thruster 

2 

design formulated for argon operation, the anode area is 0.4 m and the 
power input is 2604 W. Consequently, the steady-state anode temperature 
will be determined approximately by: 

4 1.148 x 10 llo K 

1 r 

where e is the emissivity of the anode material. Figure 9 shows the 
relationship between the steady-state temperature in °C and the emissivity. 
Since the magnets are in relatively good contact with the anode (thermally), 
the maximum allowa h1c> temperature for the anodes must be the same as for 
the permanent magnets. Figure 10 shows the thermal characteristics of 
several magnet materials, and ALNICO 8 is the obvious choice. As shown, 
550°C is the maximum allowable temperature, and the objective would be 
to limit temperature to 450°C or less to have some margin of safety 
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Figure 9. Steady-state temperature attained by radiation 
cooling as a function of surface emissivity 
(50-c;.i thruster, argon propellant). 
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REVERSIBLE TEMPERATURE EFFECTS 



ALNICO 

840C 

INDOX 

460 C 

CUNIFE 

410C 

INCOR 16 

700C 


ALNICO 

550 C 

INDOX 

1080C 

CUNIFE 

300 C 

INCOR 16 

300 C 


IRREVERSIBLE TEMPERATURE EFFECTS 


MATERIAL 

PERCENT LOSS 

-60C 

-20 C 

200 C 

ALNICO 5 

2.5 

1.0 

0.8 

ALNICO 8 

0.1 

0.1 

0.2 

INDOX 

8 

4 

0 

CUNIFE 

2.0 

1.0 

15.0 

INCOR 16 

0 

0 

2.0 


MAXIMUM PRACTICAL OPERATING 
TEMPERATURES IN AIR 


ALNICO 

540C 

INDOX 

300 C 

CUNIFE 

250 C 

INCOR 16 

250 C 


Figure 10. Magnet properties that limit thruster operating temperature 
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and also to minimize degradation. Figure 9 shows that an eraissivity 
of 0.4 or greater is required. As seen in Table 5, the emissivities 
available are on the order of 0.2 for typical materials that could be 
used to fabricate thruster anodes. Consequently, a thermal-control coating 
will be necessary to achieve satisfactory temperatures. Two thermal- 
control coatings are available that provide a minimum emissivity of 0.8 
and have been qualified for space applications. The properties of these 
thermal-control coatings are described below: 

Thermal Control Finishes (From HAC "Process Manual") 

Specification: HP 4-190 

I 

Application: This flat white coating is a high- temperature 

650° C resistant, passive thermal- control coating suitable for 
use in a space environment. The coating provides low solar 
absorptance and high IR emittance. It can be used on surfaces 
that are rapidly heated to and exposed to continuous service 
at 650°C. However, volatile material will be released at 
the latter temperature. The coating has moderate resistance 
to UV degradation (i.e., change in solar absorptance)'. 

Because of limited testing and 425°C curing temperature , 
current use of this coating is restricted to nickel-plated 
metals and CRES steels. 

Material: HMS 14-1789 

Typical Callout: Paint Per HP 4-190 


Specification: HP 4-189 

Applications: This flat black , high-temperature (650°C) 

resistant , passive thermal control coating is suitable for 
use in a space environment. The coating provides a high 
solar absorptance and high IR emittance. Because of limited 
testing and 425°C curing temperature, current use of this 
coating is restricted to nickel-plated and CRES steels. 

Material: HMS 15-1766 

Typical Callout: Paint Per HP 4-189 


By using one of these finishes on the external anode and pole-piece 
surfaces, steady-state temperatures on the order of 350 to 400°C should 
be achievable. 
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Table 5. Properties of Potential Thruster 
Fabrication Materials 


Thruster Part 
Candidate Material 

HE9 

Melting 
Point, °C 

C N 

a, “C" 1 

Anode 

Aluminum 

2.22 

646 

0.15 

23 x 10“ 6 

Copper 

3.91 

1083 

0.15 

17.6 x 10~ 6 

Beryllium 

1.51 

1285 

0.10 

11.5 x 10“ 6 

Magnesium 

0. 76 

566 

0.10 

28.8 x 10' 6 

Austenitic CRES 

0.16 

• 1400 

0.27 

18.4 x 10 -6 

Pole 

1010/1020 steel 

0.47 

1500 

0.27 

13.5 x 10~ 6 

Electrode 

Molybdenum 

1.46 

2610 

0.15 

4.9 x ,0 -6 

Titanium 

0.07 

1650 


9.2 x 1C" 6 

Magnet 

Alnico 




13.0 x 10 -6 


where: k is thermal conductivity (in W/cm°C) 


is normal total emissivity (dimensionless). 

N -1 
a is thermal expansion (°C ). 


7174 






The conceptual approach for achieving heat rejection and tractable 
equilibrium temperatures forms the basis for examining the mechanical/ 
structural design of the thruster. At the outset, noncircular, oblong 
thruster cross sections were investigated to achieve larger beam areas 
without exceeding the limitations of available molybdenum sheet material 
(necessary to fabricate the ion-optics grids). The ion-optics sub- 
assembly is unquestionably the most critical subassembly with regard 
to dimensional stability under variations in temperature. The inter- 
electrode spacing must be maintained at about 0.75 mm ± 0.2 mm for all 
temperatures within the thruster operating range. All structural and 
mechanical design features must be tailored to enable this type of 
dimensional stability in the ion-optics subassembly. Consequently, the 
ion-optics subassembly was analyzed extensively (by finite element 
analysis) to determine which factors have the most effect on inter- 
electrode spacing. 

Structurally, the 30-cm ion-optics and the 50-cm ion-optics are 
nearly identical. The 50-cm ion-optics was modeled by increasing the 
radial dimensions of the 30-cm ion-optics EASE-2 structural analysis 
program. The axial dimensions, and thus the dish depth, were not changed. 
This resulted in a proportionately shallower dish electrode structure. 

The radius of curvature of the spherically dished portion was increased 
from 0.5 m (20 in.) to 1.3 m (51 in.). Both two- and three-grid designs 
were modeled by this means. Isometric views of the models are shown in 
Figures 11 and 12. 

The thin (0.4 mm) grids of both 30-cm and 50-cm ion-optics subas- 
semblies are attached to thicker (1.9 and 2.5 mm) support rings around 
their periphery. The screen grid support ring (2.5 mm thick) is then 
attached to the mounting ring. This mounting ring consists of a thick 
(5.1 mm) ring forming the base with a thin wall (0.1 mm) short cylin- 
drical section at the inside diameter. The cylindrical section is slotted 
about two-thirds the distance from the open end towards the base on both 
sides of 12 attachment points for joining the screen grid support ring 
to the mounting ring. The slots, located 0.6 cm from the attachment 
points, form a flat spring column 1.2 cm wide and 2.3 cm high which may 


44 



ORIGINAL PAGE 13 
0F POOR QUALITY 


* 67 3 - 44 



F1 re 11. Isometric view of the computer modeling for 
8 a 2-grid 50-cm thruster ion-optics assembly. 


V673-46 



Figure U. Isometric view of the computer modeling for 
a 3-grid 50-cm thruster ion-optics assembly. 
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flex radially and is stiff tangentially through its width. The accel 
(and for three-grid, the decel) are attached to the base of the mounting 
ring. They are electrically isolated by alumina insulator assemblies 
mounted to the base with flat spring assemblies, similar in form and 
fvn -tion to the flat spring columns in the mounting ring, between the 
insulators and the grid support rings. 

The support rings and mounting ring are modeled as solid elements 
and the insulator as pipe elements; the flat spring flex member is modeled 
to support axial loads only. For the analytic model, the grids are 
attached to the mounting ring with rigid beams (at the base only), and 
the slots in the mounting ring next to the shoulder are included in the 
model. The modeling is three dimensional. A 30° section, shown in 
Figure 11, was modeled as defined by symmetry of the structure. 

The center nodes of the accel and the screen electrodes are allowed 
to move only in the axial direction. The bottom of the mounting ring 
is considered to be placed on rollers. 

The grid and support rings are molybdenum, the mounting ring is 
titanium alloy, the insulator is alumina ceramic, and the flex member is 
stainless steel. 

Thermal loading was varied to correspond to temperatures that were 
measured for operation of the 30-cm thruster at beam currents of 2 A 
and 6 A, and for a temperature 100°C greater than those for operation at 
a beam current of 6 A, shown in Figures 13 through 15. In modeling the 
two-grid ion-optics assembly, temperature data were obtained from 
NASA LeRC tests for thruster operation at 2 A: temperature data for 

thruster operation at 6 A were obtained from Hughes Advanced Technology 
Program. Since the data originated from tests of two-grid ion-optics 
assemblies, temperature data for three-grid ion-optics assemblies were 
derived by assuming the accel temperature of the two-grid assembly to be 
the decel temperature of the thvee-grid assembly; accel temperature of 
the three-grid assembly was then assumed to be higher than the screen 
grid temperature by one-half the differential between the temperatures 
of the screen and decel grids. 
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Figure 13. Cross section of the ion accelerator grid assembly for the 

J-series 30-cm thruster showing temperatures used in stress 

analysis J = 2 A, 
b 

8996- 17 R1 



Figure 14. Cross section of the ion accelerator grid assembly showing 
temperatures representative of 1.4 W/cm 2 thermal input 

J b = 6 A ' 
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Figure 15- Ion-optics model showing temperatures assumed for condition 

that operating temperature would be 100°C greater than for 

J * 6 A. 
b 

The results of these analyses show good stability of all of the 
electrodes with interelectrode spacing varying only due to thermal 
expansion. There was no apparent deformation of the electrodes from 
electrode supporting rings, indicating that the. flex members connecting 
the molybdenum electrode supporting ring to the titanium electrode 
mounting ring are decoupling the effects of differential thermal expan- 
sion between the two rings. 

The accel electrode for a 50-cm-wide by 100-cm-long oval thruster 
was modeled and analyzed. A quadrant of the accel electrode was modeled 
using the two axes of symmetry as boundaries. Motion at the intersection 
was limited to axial direction only. The bottom of the mounting ring 
was mounted on rollers, free to move in two directions but not axially. 
Flex member supports connected the support ring with the mounting ring. 

A plane view of the accel electrode model is shown in Figure 16. The 
temperatures were the same as those for the round electrodes with the 
temperature at the center of the semicircular end assumed to be constant 
along the center of the rectangular part of the oval electrode. 
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Figure lb. Plan view M)-cm x 100- cm ova 1 aocel electrode 
mode 1 . 


The result of the analysis for load ease 2 (;K)-om olootrodo 
temperatures at b A beam) is shown super imposed on the room-temperature 
model in Figure 17. The room-t emperat ure model is shown in this figuta 
as broken lines. Thermal displaeement of the nodes in the rectangular 
part ot the grid is nonunitorm, with the hotter central axis expanding 
more than the cooler outside edges. This causes the circular quadrant 
to ho distorted such that tin* point of maximum displacement parallel 
to t he long axis ot the electrode is along a radial at 10" from the 
center ot the circular quardrant , not along the long axis as would he 
expected. i'onsequent 1 v, t he symmetrical expansion necessary to preserve 
electrode .".pert ure alignment cannot be assured. 

The quasi-radial lines hovond tin' outline ot the grid nodes are 
the two ends ot the tlex members. The outer -most end is the end that is 
attached to the insulators. The inner ends are attached to the accel 
stiffening ring. The flex members along the two axes are the least 
attected by d i I t erenl ia 1 thermal expansion between the titanium electrodi 
and molybdenum stiffening rings. The other tour flex member supports 
would have to be able to tlex in two directions and carry axial loads 
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only. Consequently, non-circular thruster cross sections are not con' 
sidered to be tractable from an ion optics assembly point of view (at 
least for large thrusters that are operated at high temperature). 
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Figure 17. Displacements of model elements in thermal expansion. 
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SECTION 3 


RECOMMENDED CONCEPTUAL DESIGN 

The analyses and computations described in the previous section lead 
to the recommendation of a 50-cm— diameter multipole thruster configured 
as shown in the layout drawing of Figure 18. The main design specifica- 
tions and performance characteristics are listed in Table 6. The design 
features that distinguish this conceptual design from the NASA/Hughes 
30-cm mercury ion thruster are listed as follows: 

• Modular, multipole magnetic confinement discharge chamber 

• Three-electrode i^n-optics assembly for high-current, 

low-1 operation 
sp ^ 

• Anode potential support structure to facilitate heat 
rejection 

• Multiple cathodes to build on state-of-the-art 
hollow-cathode technology 

• Plenum-type distributed electron injection with 
electromagnets for control 

• Major subassemblies electrically isolated to enable 
choice of bias 

• Propellant flow controlled with piezoelectric valves 

• Thruster assembled from major subassemblies. 

The major structural element in this design is the angle cross- 
section ring at the intersection of the end and the side walls. This 
ring is fabricated of austenitic stainless steel for its non-magnetic 
characteristic. The angle cross section supplies the stabili:y required. 
Not only is it the main structural element, but it also serves as the 
main propellant plenum and is shown with four main isolators directly 
feeding into the upstream wall of the ring. The arms of the angle are 
partially closed by the external magnetic poles of the end and side 
wall. Two annular slots are provided by the inclusion of the internal 
corner magnetic pole, which also has an angle cross section. Propellant 
distribution uniformity will be provided by a fine wire mesh cloth ring 
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Figure 18. Layout drawing of the proposed thruster diagram. 
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Tabic 6. Thruster Design Specifications and Predicted 
Performance Characteristics 


specification or Character! st 1c 

Design 

Value 


Argon Propellant 

Xenon Propellant 

Specific impulse, sec 

6,076 

3,530 

Thrust, N 

0.5 

0.5 

Beam diameter, cm 

50 

50 

Discharge chamber length, cm 

16.2 

16.2 

Beam voltage, V 

932 

716 

Discharge voltage, V 

45 

34 

emission current, \ 

86.5 

63.1 

Beam current, A 

18.9 

l5.0 

Total power input, W 

21,500 

11,453 

Electrical efficiency, % 

81. 9 

81. 3 

Propellant efficiency, uncorrected, % 

93. 5 

122.6 

Thrust correction factor 

0.947 

0.869 

Thruster efficiency, % 

68.6 

75.2 

Screen grid wear out, hr 

6,130* 

34,000 

Specific discharge power, eV/ion 

206 

165 

Total extraction voltage, V 

1,550 

1,795 

Ratio of beam voltage to total 
voltage 

0.6 

0.4 

Number of cathodes 

4 

4 

Total propellant flow. A, equivalent 

20. 2 

10.6 

Cathode propellant flow, each. A, 
equivalent 


1.4 

Neutralizer propellant flow. A, 
equivalent 

1.5 

1.0 

Propellant efficiency including 
neutralizer, % 

8/.0 

112.1 

Electrical efficiency including 
neutralizer, % 

80.6 

79.6 

Thruster efficiency including 
neutralizer, % 

62.9 

67. 3 

‘C.n be extended to 24,500 by adding 0.26% (mass) nitrogen to the argon 
propel lant ,15,1b 
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suspended between the two external magnetic poles aided by baffles 
which will deflect the incoming propellant stream from the isolators. 

The end wall anodes will be attached to the external magnetic pole. 

The end wall anode will be machined in two rings, separated bv the 
electron plenum, an annular chamber containing the multiple hollow 
cathode. As shown in Figure 18, the electron plenum is electrically 
isolated from the anode and cathodes so that the electron plenum can be 
maintained at keeper potential. For maximum flexibility, the prototype 
thruster would be constructed so that the cathode, cathode keeper, and 
the electron plenum can all be isolated from each otherjand their 
potentials can be changed during tests without breaking chamber vacuum 
or touching the thruster. 

The end wall magnetic poles can be easily formed from strips of 
1010 to 1020 steel. The external poles will be attached directly to 
the anode, while the internal poles will be electrically isolated with 
ceramic insulators. A maximum gap betweei matching external and 
internal poles of 0.12 cm will be maintained by locating the external 
poles in annular slots cut in the upstream surface of the end wall 
anodes. The poles will be segmented to allow for differential thermal 
expansion between the copper or stainless-steel anode and the steel poles. 

The side wall is modular, consisting ot alternating flat steel 
magnetic pole rings and the cross section anode rings. This construc- 
tion provides the flexibility to change the chamber length in increments 
of 2.7 cm, which is the nominal cyclic length, for the multipole thruster 
as developed bv the researchers at CSU. Although the analysis described 
earlier shows that xenon and argon will require different chamber lengths 
to obtain maximum efficiency, a single discharge chamber length of 16.2 cm 
Is recommended. Chamber length can he varied by adding or removing poles 
and anodes as determined bv testing. 

The electrode assemble Is supported by sheet metal columns attached 
to the outer-ring propellant plenum. An intermediate ring with a chan- 
nel cross section is also shown as a support for the insulation isolating 
the electrode assembly from the thruster. This intermediate ring and 
t lie electrode mounting ring will be fabricated from titanium alloy to 
minimize the possibility of distorting the electrode mounting ring 


due to differential thermal expansion. The electrode mounting ring 
is a flat ring of approximately 0.63 cm. thickness to present the least 
profile to the anode. This cross section may not provide sufficient 
rigidity or stability, which would require the addition of a rim or 
rims to increase the section modulus. As in the design of the 30-cm 
mercury thruster, the thruster is composed of major subassemblies for 
ease in fabrication and to enable refurbishment. The major subassemblies 
are listed as follows: 

• Three-electrode ion optics and support ring 

• Outer multipole discharge-chamber boundary and magnets 
(separate designs for xenon and argon) 

• Rear wall multipole discharge-chamber boundary and 
magnets 

• Electron injection plenum with hollow cathodes 

• Propellant injection plenum 

• Propellant valve and isolator 

• Ground screen 

• Neutralizer. 


A. POWER PROCESSOR REQUIREMENTS 

The power processor requirements are shown conceptually in the 
schematic shown in Figure 19 (note that the power supplies have dual 
outputs). There are four major power source requirements — propellant 
control, discharge, screen/accel, and neutralizer. A power supply has 
not been included for the electromagnets used to adjust the magnetic 
field in the electron injection gap because it is thought that permanent 
magnets will replace the electromagnet after the initial thruster 
development. The control concept for controlling discharge voltage, V^, 
and for balancing cathode emission is shown in Figure 20. The gas flow 
through the valve supplying each hollow cathode is controlled in 
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Figure 19. Block diagram showing proposed thruster power processor 
requirements. 
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Figure 20. Block diagram showing four-cathode control concept, 
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proportion to the output of the control circuit, <j>^, according to the 
expression 


*1 " (V D - V REF + V C ‘ ->l C i )G t • 

where the symbols are defined in Figure 20. The power input to the con- 
trol circuitry for the piezoelectric valves should be negligible with 
respect to the other power inputs and has not been included in the effi- 
ciencies quoted previously. The nominal power supply requirements for 
the specifications listed earlier in Table 6 are shown in Table 7. It 
shoold be recognized that the voltage and currents listed in Table 7 are 
not the recommended maximum ratings required for the thruster power sup- 
plies. Power supply specifications would have to be determined on the 
basis of reserve margins that are considered adequate to the applications 
of the thrust systems in addition to the nominal power levels considered 
here. 


B. TECHNOLOGY AREAS THAT REQUIRE VERIFICATION 

The proposed design as described in the previous section is based 
on systematic extrapolation of analytic models for describing multipole 
thruster operation. Many of these concepts have been correlated with 
experimental data; however, the proposed design includes several concepts 
that have not yet been demonstrated. There are several technology areas 
that apply to large multipole thrusters in general, and some that are 
more relevant for use with a specific propellant. Some of the general 
technology areas that require verification are 

• Verification of the uniform discharge plasma conditions that 
are assumed in the performance analysis 

• Demonstration of discharge chamber operation with two or 
more hollow cathodes 

• Balancing of emission between hollow cathodes by control 
of gas flow 

• Determination of operating criteria for long-life operation 
of hollow cathodes on inert gases 
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• Verification of dimensional stability of ion optics grids 
that have longer radius of curvature (than 30-cm thruster 
grids) 

• Verification of dimensional stability at higher operating 
temperatures (than demonstrated for 30-cm thrusters) 

t Demonstration of the effectiveness of thermal control 
coatings for reducing the operating temperature of an 
ion thruster 

• Verification that the value of /Bdx used enables 
achievement of the discharge loss parao"*te’* assumed in 
the analysis 

• Verification of the magnet requirement specification for 
pole-piece geometries with an "air gap" 

• Development of a valid thermal model. 

Many of the items in the list above are not really dependent on thruster 
size or the propellant used (e.g. , multipole cathode operation, grids 

Table 7. Power Supply Requirements for the Proposed 
50-cm-Diameter Inert Gas Ion Thruster (s) 


Power Supply 

Nominal Requirements 

Voltage, V 

Current , A 

Power , W 

Argon propellant 




Screen 

900 

19.0 

17,100 

Accel 

700 

0.05 

35 

Discharge 

45 

105.4 

4,743 

Neutralizer 

15 

4 

60 

Xenon propellant 




Screen 

700 

13 

9,100 

Accel 

1,100 

0.03 

33 

Discharge 

34 

76 

2,584 

Neutralizer 

15 

*> 

30 
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with increased radius of curvature) and are being investigated under 

i 12,13,14 

several other programs. 

For operation of the proposed thruster design on xenon propellant, 
several unique technology issues are raised: 

• Feasibility of simultaneously achieving the relatively high 
propellant utilization efficiency predicted theoretically 
and the relatively low discharge voltage required for 
limiting multiple ionization 

• Dimensional stability of three-grid ion optics assembly 
(operated at R ■ 0.4) 

Operation of the thruster designed for use with argon propellant raises 
still other technology issues 

• Achievement of adequate lifetime with the high discharge 
voltage required (45 V) 

• Feasibility of reducing discharge chamber sputtering 
by adding nitrogen to the propellant 

• Effect of nitrogen propellant additive on hollow 
cathode operation and lifetime 

• Feasibility of achieving long operating lifetime with 
hollow cathodes operated at greater than 20-A emission 
current . 

Experiments will be required to settle these issues. 


SECTION 4 


CONCLUSIONS 

The design analyses performed during this study have shown that an 
ion thruster of ^50 cm in diameter will be required to produce a thrust 
of 0.5 N using xenon or argon as propellants, and operating the thruster 
at a specific impulse of 3530 sec or 6076 sec (respectively). A multi- 
pole magnetic confinement discharge chamber was specified because tne 
thruster can (as has been demonstrated experimentally) be scaled to 
other values of thrust or beam diameter without the requirement for 
appreciable empirical optimization. Because of this fact, an analytic 
description of the thruster operating characteristics is considered valid. 
Achievement of the thruster efficiency goals of 68% for argon and 74% for 
xenon was determined to be the most demanding of the thruster performance 
goals (limited by basic physical lawn). Achieving this goal requires that 
the value of the surf ace- to-volume ratio of the primary electron region 
of the discharge chamber be adjusted to obtain an acceptable balance 
between the ion production cost (power) and the resultant propellant 
efficiency. 

Other design challenges were identified as listed below: 

• Heat rejection at high power density 

• High ion beam current at relatively low net 
acceleration voltage 

• Long operating lifetime at high discharge plasma 
density (erosion by ion sputtering) 

• Dimensional stability under severe thermal load 

Engineering innovations have been deviled to answer each of these 
challenges. The heat generated can be radiated by providing good thermal 
paths through the thruster body to high emissivity surfaces (treated with 
thermal control coatings). High ion beam currents can be generated at 
low net acceleration voltages (high perv"ani_el by using three-grid ion 
optics. Long operating lifetime can be achieved by operating at low 
discharge voltage and/or adding nitrogen to the propellant to reduce 
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sputtering rates. Dimensional stability can be achieved by using flexible 
members in appropriate places to permit predictable elastic deformation 
of thruster components without appreciably changing critical relative 
dimensions (such as interelectrode spacing). A circular thruster (ion 
beam) cross section was found to be necessary for ensuring predictable 
elastic deformation under severe thermal loads. 

The conceptual design that was formulated includes several design 
features that require experimental verification. The more important of 
the untested design approaches are listed below: 

• Achievement of dimensional stability of dished ion 
optics grids that have a 1.3-m radius of curvature 

• Provision of a satisfactory magne :ic field distribution 
in the anode/pole-piece region when an "air gap" is 
introduced into the pole-piece magnetic circuit 

• Achievement of uniform discharge plasma conditions by 
operating the thruster with multiple cathodes (balancing 
current between cathodes) 

• Feasibility of reducing the discharge chamber sputtering 
by adding nitrogen to the propellant 

• Feasibility of achieving enhanced thermal rejection through 
use of thermal control coatings. 

The design features listed above are considered to be reasonable 
extrapolations on exploratory work in progress and, therefore, are not 
great developmental risks. Assuming that the engineering details for 
these design innovations can be successfully implemented, the resultant 
thruster will fulfill the design goals set forth under this contract 
(exclusive of the neutralizer). Using state-of-the-art technology for 
hollow cathodes operated on inert gases to predict neutralizer propellant 
requirements leads a severe performance penalty; however, an extrapo- 
lation of achievable performance improvements would be pure conjecture at 
this time. Consequently, detailed design specifications for the neutral- 
izers have been left as an open item under this study. Similarly, 
propellant electrical isolators and gas control were considered in 
concept only . 
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At the outset, and early in the study, the contractual goals seemed 
unattainable within reasonable extrapolations on proven state-of-the-art 
technology and analyses. It now appears that the conceptual thruster 
design that has evolved during this study can be developed to meet or 
exceed these goals without the requirement for any major technological 
" breakthroughs . " 


63 



REFERENCES 


1. H.R. Kaufman, "Experimental Investigation of Argon and Xenon Ion 
Sources*" Annual Report, CR-134845, Colorado State University, 

June 1975. 

2. H.R. Kaufman, "Inert Gas Thrusters," Annual Report, CR-135100, 
Colorado State University, July 1976. 

3. L.A. Rehn, "Argon Hollow Cathode," NASA CR-135102. Thesis for 
Grant NSG-3011, November 1976. 

4. H.R. Kaufman, "Industrial Ion Source Technology," NSG-3086, 

NASA CR-135149, November 1976. 

5. H.R. Kaufman, "Inert Gas Thrusters," NSG 3011, CR-135226, Colorado 
State University, July 1977. 

6. G.C. Isaacson, "Multipole Gas Thruster Design," NASA CR-135101, 
for Grant NSG-3011, June 1977. 

7. H.R. Kaufman and P.J. Wilbur, "Scaling of Mercury and Gaseous 
Propellant Ion Thrusters to Large Sizes," AIAA 78-667. Presented 
at AIAA 13th International Electric Propellant Conference, San 
Diego, California, April 1978. 

8. L. Rehn and H.R. Kaufman, "Correlation of Inert Gas Hollow Cathode 
Performance," AIAA 78-707, Colorado State University. Presented at 
AIAA 13th Electric Propellant Conference, San Diego, California, 
April 1978. 

9. H.R. Kaufman, "Industrial Ion Source Technology," NSG-3086, 

CR-135353, November 1977. 

10. P.T. Wilbur and H.R. Kaufman, "Scaling Relationships for Mercury 

and Gaseous Propellant Ion Thrusters," AIAA Paper 78-667, April 1978. 

11. H.R. Kaufman, "Inert Gas Thrusters," NSG 3011, NASA CR-159527, 
Colorado State University, 1978. 

12. H.R. Kaufman and R.S. Robinson, "Plasma Processes in Inert Gas 
Thrusters," Paper 79-2055, Princeton University, AIAA/DGLR 14th 
Electric Propulsion Conference, November 1979. 

13. G.R. Longhurst and P.J. Wilbur, "Plasma Property and Performance 
Prediction for Mercury Ion Thrusters," Paper No. 79-2054, presented 
at Princeton/AIAA/DGLR 14th Electric Prop-tlsion Conference, 

Princeton, NJ, November 1979. 


65 


i •lu.wtJfcNia PAGE BLANK NOT FILMED 



14 . 


J.R. Bropl.y and P.J. Wilbur, "Electron Diffusion through the Baffle 
Aperture of a Hollow-Cathode Thruster," Paper No. 79-2060, presented 
at Prineeton/AIAA/DGLR 14th Electric Propulsion Conference, Prim "ton, 
NJ, November 1979. 

15. P.J. Wilbur, "A Model for Nitrogen Chemisorption in Ion Thrusters," 
Paper No. 79-2062, presented at Princeton/AIAA/DGLR 14th Electric 
Propulsion Conference, Princeton, NJ, November 1979. 

16. M.A. Mantenieks and V.K. Rawlln, "Sputtering in Mercury Ion Thrusters," 
Paper No. 79-2061, presented at Princeton/AIAA/DGLR 14th Electric 
Propulsion Conference, Princeton, NJ, November 1979. 



APPENDIX A 


THRUSTER PERFORMANCE CHARACTERIZATION 


Optimum discharge-chamber performance is normally obtained near 
the "knee" of the discharge-chamber performance curve. This is because 
operation at higher discharge losses will substantially increase the 
power loss for a much smaller increase in utilization. Conversely, 
operation at lower discharge losses will decrease propellant utiliza- 
tion more than the power loss will be decreased. The exact point near 
the knee that will prove optimum for a given set of thruster operating 
specifications will depend on the required operating conditions. To 
simplify the iteration procedure, however, a single knee point will be 
used for this calculation. The error associated with this assumption 
should be less than the uncertainty in the discharge-chamber performance, 
when the latter is projected from a limited number of tests. For ease 
of calculation and rapid hardware development, the multipole type of 
discharge chamber was assumed. 

A-l 

Using earlier correlations obtained with argon and xenon, 
together with the screen open-area fraction of 0.65 for the ion optics 
assumed, the knee discharge losses (eV/ion) used here are 


E = 73 (Ap/^) (A l) 

E = 67 (A p /A b ) (Xe) . 
The corresponding neutral losses (A-equiv.) are 


(A. la) 
(A. lb) 


J (1 - n.) * 3.4 A /(? /A ) (Ar) 
o i o p p 


(A. 2a) 


J (1 - n.) = 0.41 A /(y /A)(Xe) , 

o i o p p 


(A. 2b) 


where A is the effective open area (<}) A, ) of the ion optics (m ), and 

O 3D 

yp and Ap are the volume and outside area of the primary electron region 
(m^ and m^) . The effective open area, both her and in ion optics 
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designs with accel holes significantly smaller than the screen holes, 
is the accelerator open area times the Clausing factor for the length- 
to-diameter ratio of the accelerator holes used. Since a simple cylindri- 
cal design was assumed for the discharge chamber, the primary electron 
region was assumed to be the cylinder that would just fit within the 
anodes, pole pieces, and screen grid. A simpler expression was suffi- 
cient for the correlation because only a single chamber diameter and 
optics design were used. 

A convenient starting point for the ion optics design is the 
screen hole diameter. In general, the smaller the screen holes, the 
larger the current capacity of the ion optics. Several experimental 
studies have shown, though, that screen holes smaller than about 2 ram 
in diameter depart substantially from expected space-charge-flow per- 
formance. A diameter of 1.9 mm is commonly used and is assumed here 
to be the minimum diameter. 

The choice of accelerator hole size is a compromise between the 
current-carrying capacity and neutral loss, both of which increase as 
hole diameter is increased. To a first approximation, these varia- 
tions are not dependent on the use of a two- or three-grid ion optics 
design or on the value of net-to-total voltage ratio, R, used. The 
beam current capacity is proportional to the normalized perveance (N.P.), 
while the neutral loss is proportional to the effective open area of 
the accelerator grid. A performance index can therefore be obtained by 
dividing the normalized perveance by the effective open area of a single 
aperture. The effective open area also includes a Clausing factor, so 
a typical ratio of accel thickness to screen hole diameter of 0.2 is 

assumed. Expressing the accel hole area as a fraction of screen hole 

2 

area, the performance index is N.P./(d /d ) K . Using two-grid experi- 

a s c 

mental data for an R of 0.7, a plot of this index against accei-to-screen 

hole diameter ratio, d /d , is shown in Figure A-l. The normalized per- 

a s 

veance used was for mercury propellant, but the d /d value for the maxi- 

a s 

mum would be the same for other propellants. The significant point 
about the data shown in Figure A-l is the maximum near the data point 
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for d Id ■ 0.64. More detailed ion optics data would probably give a 

3 8 

slightly different value. In addition, the exact d /d would shift 

as 

slightly for different operating conditions. For the purpose" cf this 
study, though, a value of d /d * 0.64 should be sufficiently precise. 

8 S 

For the maximum span-to-gap ratio, hence maximum thruster diameter, 
some form of ion optics using dished molybdenum grids will be assumed. 

As described in a preliminary study of large thrusters, a maximum span- 
to-gap ratio of about 600:1 can be assumed for this approach. This 
maximum span-to-gap ratio is essentially determined by thermal/mechanical 
effects and will undoubtedly vary somewhat with size, power level, and 
the specific mechanical design used. As with other assumptions, however, 
the use of a single limit value is appropriate for this analysis. 

Another ion-optic parameter that must be selected is grid spacing 


relative to screen hole diameter, £ / d . We are not just interested in 
physically large thrusters: we also want them to have high performance. 

For high current and power densities, it is necessary that the grid 
spacing be small. As shown in Figure A-2, the benefits of small £ 's 

O 

are limited for £ /d g <<: 1* The best simple parameter for acceleration 

distance between a screen grid and an accelerator grid is the effective 

acceleration distance, £ , defined by the equation given in the sketch. 

The ratio £ / d , decreases with decreasing £ / d , but the decrease becomes 
e s g s 

very small as £ /d approaches zero. A tradeoff is indicated because a 
small ^ e /d g is desirable for high current density, but the tolerance 
and electrical breakdown problems become critical as £^/d g approaches 


zero. 


* 

9573-5 



d a /d s 


Figure A-l. Ion optics performance index as a func- 
tion of d /d . 


a s 
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Figure A-2. 

Normalized acceleration 
length as a function of nor- 
malized interelectrode 
spacing. 


As a partially arbitrary assumption, a value of 0.5 was selected 

for H / d . The argument given above for a small l /d is most persuas- 
8 s 8 s 
ive for low specific impulses, where the maximum electric fields are not 

a limit, even though the smallest permissible screen hole diameter is 
used. There is, though, a thermal expansion argument that is valid at 
all specific impulses. For the maximum span-to-gap ratio, thermal expan- 
sion would be a larger fraction of screen-hole diameter if a larger 

value of ? /d were used. Specific designs (particularly high specific 
8 s 

impulses at less than maximum span-to— gap ratio) may benefit from larger 

values of Z /d , but the value selected appears to be a reasonable choice 
8 s 

for all designs. 

-9 

For ion-beam calculations, a normalized perveance of 3 x 10 (Ar) 
appears to be a reasonable value. That is high, but not so high as to 
risk serious accelerator impingement. For xenon, the normalized per- 
veance was corrected by the square root of atomic mass. Backstreaming 
calculations for the assumed geometry indicate a maximum R of 0.8, 
while the minimum R was assumed for two grids to be 0.7 (for 15° half 
angle) and for three grids to be 0.17 (same 15° half angle, i * £ d ). 
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For three grids, the deceierator grid should also have holes about 0.83 
times the screen hole diameter. The deceierator grid spacing was assumed 
equal to the accelerator grid spacing because, for the largest thrusters, 
both spacings should be subject to roughly the same limit in span-to-gap 
ratio. 

For minimum screen hole diameter, then, the minimum gap is 0,5 x 
1.9, or 0.95 nan. A safe limit on maximum electric field appears to be 
2000 V/mm, although further breakdown data with operating thrusters may 
show that this is too small for the smaller gaps. This 2000 V/ram limit 
was used for all calculations made here. Then for the smallest grid 
gap, a total voltage of 190U V (0.95 mm x 2000 V/mm) is permissible. 

The gap can be increased to above the 0.95 mm value if required to (1) mai 
tain the permissible span-to-gap ratio of 600:1 or (2) reduce the scieen- 
accel electric field to 2000 V/mm. 

For the preliminary calculations shown here, neither an off-axis 
thrust loss nor a neutralizer propellant loss was included, and the only 
power loss was assumed due to the discharge. 

To avoid an unnecessary iteration, a performance calculation at a 
given specific impulse was started by assuming a propellant utilization. 
Tiie net voltage necessary to obtain the required average propellant 
velocity was then calculated; then the complete ion optics performance 
consistent with span-to-gap and electric field limits was calculated. 

Next, the discharge chamber length (or depth) required for the assumed 
utilization was calculated. (The neutral loss is actually the parameter 
of interest for this length and is available from the assumed utilization 
and the calculated beam current. It should be noted that it is possible 
to initially assume a propellant utilization that is too high to be 
obtained by any discharge chamber length.) With chamber length known, 
the discharge loss and thruster efficiency were calculated. 

Double-ion production was considered in this analysis, although 
not in detail. From the information presented earlier in this report, 
it was assumed that low discharge voltages would reduce double-ion 
production to a level consistent with both long life and high perform- 
ance. The present understanding of electron diffusion is also assumed 
to be sufficient to permit reaching the necessary low discharge voltages. 


A. CALCULATED PERFORMANCE 

Calculated performance data are shown for argon propellant In 

Figures A-3 through A-10 and for xenon propellant In Figures A-ll through 

A-16. These data are also shown in Tables A-l through A-14. 

In Figures A-3 and A-10, performance is shown over a range of 

specific impulse for a thruster with a 0.5-m beam diameter. At each 

specific impulse, the propellant utilization was iterated to the nearest 

0.001 to maximize thruster efficiency. The beam diameter of 0.5 m is 

large, but not so large that the minimum grid spacing is Increased 

above 0.95 mm due to span-to-gap ratio (d, /£ ). The maximum efficiency 

S 

also required that the ion optics for each specific Impulse be optimized, 

although this optimization did not require any iteration. 

Optimized in this manner, the data in Figures A-3 and A-10 show 

the minimum net-to-total voltage ratio, R, at low specific impulses, the 

maximum of 0.8 at high specific impulses, and intermediate values in 

between. At low specific impulses (below about 2000 sec for Ar and 

below about 3000 sec for Xe), the minimum spacing of 0.95 mm is used, 

but the screen-accel electric field is less than the maximum of 2000 V/mm. 

The use of the minimum permissible R is necessary to obtain maximum 

beam current which, in view of neutral-loss theory, is required to 

obtain maximum utilization and thruster efficiencies. 

In the intermediate specific impulse range (about 3000 to 7770 sec 

for Ar and about 2000 to 4620 sec for Xe), the maximum beam current is 

obtained by keeping the minimum gap dictated by the span-to-gap ratio 

and decreasing the total voltage to agree with the electric field limit. 

This conclusion can be shown by Child's law, with current varying as 
3/2 2 

AV /£ , Using a gap larger than the minimum while maintaining maxi- 

t e 2 3/2 

mum electric field will result in increasing more than AV fc , thus 

decreasing maximum current. With the gap fixed at the minimum and the 

electric field at the maximum, the extracted ion current is also a 

constant in this specific impulse range (see Figures A-3(b) and A-10(b)). 

In the high specific impulse range (above abouL 7770 sec for Ar and 

about 4620 sec for Xe), the net voltage is so high that R will incre*se 
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e) DISCHARGE CHAMBER LENGTH 



f) THRUSTER EFFICIENCY 
Figure A-3. Continued. 
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c) THRUSTER EFFICIENCY 


Figure A-4. Effect of propellant utilization on perform- 
ance. Propellant, argon; beam diameter, 

0.5 m; specific impulse, 2000 sec. 
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j) THRUST 



0 5 0.6 0.7 0.8 0.9 1.0 

PROPELLANT UTILIZATION. n u 

c) THRUSTER EFFICIENCY 

Figure A-5. Effect of propellant utilization on per- 
formance. Propellant, argon; beam diam- 
eter, 0.5 m; specific impulse, 5000 sec. 
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Figure A-6, Effect of propellant utilization on performance. 

Propellant, argon; beam diameter, 0.5 tn; specific 
impulse, 10,000 sec. 
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Figure A-7. Effect of beam diameter on performance. 

Propellant, argon; specific impulse, 2000 sec; 
propellant utilization, optimized for maximum 
thruster efficiency. 


78 





SPECIFIC THRUST 

THRUSTER EFFICIENCY (T/A^), N/m 2 CHAMBER L/D 


ORIGINAL PAGE (S 
OF POOR QUALITY 


9673— 13 



a) CHAMBER LENGTH-TO-DIAMETER RATIO 



b) THRUST-TO-BEAM AREA RATIO 
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BEAM DIAMETER, m 
c) THRUSTER EFFICIENCY 


Figure A-8. Effect of beam diameter on performance. Propel- 
lant, argon; specific impulse, 5000 sec; propel- 
lant utilization, optimized for maximum thruster 
efficiency. 
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Figure A-9. Effect of beam diameter on performance. Propellant, 
argon; specific impulse, 10,000 sec; propellant 
utilization, optimized for maximum thruster 
efficiency. 
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c) THRUST 


Figure A-10. Effect of specific impulse on voltage ratio, beam cur- 
rent, thrust, thrust to power ratio, discharge chamber 
length, and thruster efficiency. Propellant, xenon; 
beam diameter, 0.5 m; propellant utilization optimized 
for maximum thruster efficiency. 
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e) DISCHARGE CHAMBER LENGTH 



f> THRUSTER EFFICIENCY 
Figure A- 10. Continued. 
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Figure A-ll. Effect of propellant utilization on perform- 
ance. Propellant, xenon; beam diameter, 0.5 m, 
specific impulse, 1500 sec. 
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c) THRUSTER EFFICIENCY 

Figure A-12. Effect of propellant utilization on performance. 

Propellant, xenon; beam diameter, 0.5 m; specific 
impulse, 3000 sec. 
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b) THRUST 
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PROPELLANT UTILIZATION (t| u ). 
c) THRUSTER EFFICIENCY 

Figure A-13. Effect of propellant utilization on performance. 

Propellant, xenon; beam diameter, 0.5 ra; specific 
impulse, 10,000 sec. 
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Figure A-15. 


Effect ol beam diameter on performance. Pro- 
pellant, xenon; specific impulse, 3000 sec; pro- 
pellant utilization, optimized for maximum 
thruster efficiency. 
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tcAM DIAMETER, m 
c» THRUSTER EFFICIENCY 

Figu:«i A-16. Effect of beam diameter on performance. Propel- 
lant, xenon; specific impulse, 10,000 sec; pro- 
pellant utilization, optimized for maximum 
thruster efficiency. 
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Table 

— i - . a-,.. . ... — i 

A- 6 Argon Prop* 

llant. 

Opt imun 

—4. - i ■ ■ - * 

Propellant Utilization, 5000 sec. 



0.20 

0.762 

0.597 

3.93 

0.096 

3.59 

857. 

1900. 

0.4512 

0.950 

0.063 

236. 

0.30 

0.805 

0.627 

7.96 


8.08 

766. 

1900. 

0.4033 

0.950 

0.075 

219. 

0.40 

0.331 


13.3 

0.351 

14.4 

721. 

1900. 

0.3794 

0.950 

0.034 

207. 

0.30 

0.849 


20.0 

0.537 

22.4 

691. 

1900. 

0.3635 

0.950 

0.093 

200. 

0.57 

0.859 

0.665 

25.4 

0.689 

29.2 

675. 

1900. 

0.3551 

0.950 

0.099 

197. 

0.70 

0.866 

0.670 

34.0 

0.931 

39.7 

664. 

2333. 

0.2845 

1.167 

0.114 

194. 

O.SO 

0.870 

0.674 

41.2 

1.132 

48.5 

638. 

2667. 

0.2466 

1.333 

0.125 

192. 

1.00 

0.877 

0.679 

56.7 

1.569 

67.8 

647. 

3333. 

0.1942 

1.667 

0.147 

189. 

1.139 

0.878 

0.682 

63.5 

1.905 

82.4 

646. 

3797. 

0.1701 

1.898 

0.156 

186. 

1.50 

0.855 

0.672 

77.4 

2.120 

89.3 

681. 

4006. 

0.1700 

2.500 

0.204 

186. 

2.00 

0.836 

0.660 

86. 2 

2.320 

95.5 

712. 

4190. 

0.1700 

3.333 


190. 



Table 

A— 7 Argon Propellant, 

Optimum Propellant 

Utilization, 10 

,000 sec. 



0.20 

0.787 

0.697 

8.96 

0.127 

2.47 

3215. 

4019. 

0.3000 

2.009 

0.1S3 

414. 

0.50 

0.888 

0.797 

49.0 

0.796 

17.4 

2525. 

3156. 

0.3000 

1.578 

0.242 

288. 

0.879 

0.922 

0.835 

144. 

2.460 

55.9 

2342. 

2930. 

0.7994 

1.465 

0.296 

244. 

1.0 

0.92s 

0.833 

174. 

2.975 

67.8 

2327. 

3333. 

0.6932 

1.667 

0.324 

241. 

2.0 

0.940 

0.854 

476. 

8.290 

191. 

2253. 

6667. 

0.3387 

3.333 

0.532 

224. 

3.90 

0.949 

0.867 

1263. 

22.34 

522. 

2211. 

13000. 

0.1701 

6.500 

0.838 

209. 

5.0 

0.941 

0.S59 

1320. 

23.12 

536. 

22-9. 

13223. 

0.1700 

3.3:: 

1.1S6 

215. 

10. 

0.913 

0.830 

1541. 

26.09 

587. 

2339. 

14052. 

0.1700 

16.67 

3.154 

238. 

20. 

0.874 

0.793 

1920. 

31.07 

669. 

V 

15334. 

0.1700 

jj.33 

3.134 

265. 


Table 

CO 

Xenon Propellant 

,0.5m 

Beam Diameter, 0 

pticun Pr 

opellant 

Utilization. 


1000. 

0.613 

0.340 

1.54 

0.107 

4.90 

174. 

1025. 

0.1700 

0.950 

0.011 

140. 

1500. 

0.799 

0.402 

2.80 

0.137 

7.47 

231. 

1357. 

0.1700 

0.950 

0.019 

144. 

2000. 

0.900 

0,617 

5.33 

0.367 

12.4 

323. 

1900. 

0.1701 

0.950 

0.026 

148. 

2500. 

0.926 

0.700 

7.82 

0. 446 

12.4 

477. 

1900. 

0.2511 

0.950 

0.033 

154. 

3000. 

0.939 

0.758 

10.3 

0.528 

12.4 

663. 

1900. 

0.3517 

0.950 

0.048 

160. 

4000. 

0.953 

0.830 

16.4 

0. o94 

12.4 

1153. 

1900. 

0. 6069 

0.950 

0.0b8 

170. 

4620. 

0.959 

0.858 

71.0 

0.79b 

12.4 

1519. 

1900. 

0.7996 

0.950 

0.032 

180. 

5000. 

0.959 

0.370 

•> 7 * 

0.796 

11.4 

1779. 

222-. 

O.SOOO 

1.112 

0.091 

133. 

6000. 

0.960 

0.891 

26.3 

0.796 

9.54 

2557. 

3196. 

0.S000 

1.593 

0. 122 

200. 

7000. 

0.960 

0.904 

30.2 

0.796 

3.13 

3480. 

4350. 

O.SOOO 

2.175 

0.155 

217. 

8000. 

0.960 

0.912 

34.2 

0. 796 

7.16 

45-6. 

5662. 

0.8000 

2.641 

0. 195 

233. 

9000. 

0.959 

0.918 

38.3 

0.796 

6. 3b 

5765. 

7206. 

O.SOOO 

J.b03 

0.231 

258. 

100CC. 

0.953 

0. 922 

42.4 

0.79b 

5.71 

7132. 

3915. 

0.3000 

4.458 


280. 


Table A- 

9 Xenon 

Propellant, 0.5 

n Beam 

Oiar.eter, 

1500 sec 





0.50 

0. .07 

8.97 

0.496 

12.4 

589. 

1900. 

0. 3101 

0. °50 

0.003 

135. 


0.60 

0.450 

6.75 

0.-13 

12.4 

409. 

1900. 

0.2151 

0.950 

0.004 

136. 


0.65 

0. 467 

6.01 

0.381 

12.4 

3-9. 

1900. 

0.13)5 

0.950 

0.005 

137. 


0.675 

0.-74 

5.70 

0.367 

12.4 

323. 

1900. 

0.1701 

0.950 

0.00b 

137. 


0.75 

0.489 

3.63 

0.241 

9.0) 

2b2 . 

1540. 

0.1700 

0.950 

0.011 

140. 


0.799 

0. 402 

2. SO 

0.137 

7.47 

231. 

1537. 

0.1700 

0.950 

0.019 

144. 


0.35 

0 . 4 o 7 

2.21 

0.140 

6.21 

204. 

1199. 

0.1700 

0.950 

0.034 

152. 


0.90 

0.461 

1.85 

0.116 

5.23 

182 . 

1070. 

0. 1700 

0.950 

0.07) 

173. 


0.9) 

0.415 

1.31 

0.102 

4,74 

1*0. 

1002. 

0. 1700 

0.950 

0.144 

211. 


0.95 

0. 1)2 

2.07 

0.0 ’37 

4.45 

If). 

960. 

0.1700 

0.950 

0.316 

303. 


Table A- 

10 Xenon Propel!. in:, 0.5 

m Beam 

Diareter 

, 3000 sec . 




0.5J 

0.473 

24.8 

0.796 

9.94 

2356. 

294b. 

0.3000 

1.473 

0.003 

136. 


0.55 

0.514 


0.796 

10.9 

19-3. 

24 )4. 

0.80O0 

1.21? 

0 . 1)04 

136. 


0.62) 

0.572 

20.5 

0.7'’b 

12.4 

1518. 

1900. 

0.7939 

0.950 

0.0.14 

136. 


0.75 

0. 6f>2 

14. 7 

O.ftnl 

12.4 

10i\ 

1900. 

0.5512 

0.950 

0.00,8 

133. 


0.35 

0. 724 

11.9 

0. 5:1 J 

12,4 

815. 

1400. 

0.4292 

0.950 

0.016 

143. 
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Table A— 10 Xenon Propellant, 0.5 n Eeaa Diaaeeer, 5000 eee. (Continued) 


2 . 

V * 

* -£ 

s 3 


0.90 

0.743 

10.8 

0.551 

12.4 

727. 

1900. 

0.3828 

0.950 

0.023 

143. 

0.93 

0.757 

10.4 

0.553 

12.4 

631. 

1900. 

0.3585 

0.950 

0.041 

156. 

0.939 

0.758 

10.3 

0.523 

12.4 

668. 

1900. 

0.3517 

0.950 

0.043 

160. 

0.95 

0.756 

10.2 

0.522 

12.4 

653. 

1900. 

0.3436 

0.950 

0.063 

163. 

0.97 

0.732 

10. 3 

0.511 

12.4 

626. 

1909- 

0.3295 

0.950 

0.129 

203. 

0.98 

0.675 

11.0 

0.506 

12.4 

613. 

1900. 

0.3223 

0.9S0 

0.263 

277. 


Table A-ll Xenon Propellant, 0.5 m Beam Diameter, 10,000 sec. 


0.80 

0.788 

49.5 

0.796 

4.77 

10223. 

12735. 

0.3000 

6.392 

0.031 

150 

0.85 

0.836 

46.7 

0.796 

5.07 

9060. 

11325. 

0.3000 

5.662 

0.043 

'.57 

0.90 

0.831 

44.3 

0.796 

5.37 

8031. 

10102. 

0.3000 

5.051 

0.070 

172 

0.93 

0.907 

43.1 

0.796 

5.55 

7568. 

9460. 

0-8000 

4.730 

0.114 

195 

0.95 

0.920 

42.4 

0.796 

5.67 

7253. 

9066. 

0. SrOO 

4.533 

0.195 

238 

0.956 

0.922 

42.4 

0.796 

5.71 

7132. 

8915. 

0.3000 

4.458 

0.272 

280 

0.96 

0.922 

42.4 

0.796 

5.73 

7103. 

8S73. 

0. 5000 

4.439 

0.302 

296 

0.97 

0.905 

43.1 

0. 79t> 

5.79 

6957. 

8696. 

0. 3000 

4.343 

0,675 

496 


Table A— 12 Xenon Propellant, Optimum Propellant I’t ilitat ion, 1500 sec. 


0.20 

0.744 

0.477 

0.615 

0.0399 

1.48 

266. 

1565. 

0.1790 

0.950 

0.011 

149 

0.30 

0.769 

0.484 

1.19 

0.0786 

3.02 

249. 

1465. 

0.1700 

0.950 

0.014 

147 

0.40 

0.787 

0.489 

1.92 

0.127 

5.00 

233. 

1399. 

0.1700 

0.950 

0.C17 

145 

0.50 

0.799 

0.492 

2.80 

0.137 

7.47 

231. 

1357. 

0.1700 

0.950 

0.019 

344 

0.57 

0.807 

0.494 

3.48 

0.224 

9.42 

226. 

1330. 

0.1700 

0.950 

0-020 

143 

0.70 

0.795 

0.491 

3.72 

0. 243 

9.86 

233. 

1371. 

0.1700 

1.167 

0.027 

1 -4 

1.0 

0.774 

0.4S5 

4.19 

0.276 

10.7 

246. 

1446. 

0.1700 

1.667 

0.046 

146 

1.5 

0.749 

0.478 

4.85 

0. 315 

11.3 

263. 

1544. 

0.1700 

2.500 

0.0S3 

149 

2.0 

0.730 

0.472 

5.44 

0.349 

12.7 

276. 

1626. 

0.1/00 

3.233 

0.125 

151 


Table A-13 Xenon Propellant, Optimum Propellant Ccilitation, 3000 sec. 


0.20 

0.900 

0.725 

1.79 

0.0331 

1.93 

727 . 

1900. 

0.3323 

0.950 

0.031 

176. 

0.30 

0.920 

C.741 

3.85 

0.194 

4.46 

696. 

1900. 

0.3663 

0.950 

0.03? 

168. 

0.40 

0.931 

0.751 

6.63 

0.341 

7.42 

680. 

1900. 

0.3577 

0.950 

0.043 

163. 

0.50 

0.939 

0.753 

10.3 

0. 523 

12.4 

663. 

1900. 

0.3517 

0.950 

0. 0-.3 

loO . 

0.57 

0.94 3 

0.761 

1 J. 2 

0. o3J 

16.1 

662. 

1900. 

0.34S7 

0.950 

0.052 

1 58 . 

0.60 

0.944 

0.762 

14.2 

0.7)7 

17.4 

661. 

2000. 

0.3305 

1.000 

0.054 

153. 

0.70 

0* 

0. 764 

17.8 

0.927 

21.9 

653. 

2 3)3. 

0.2S23 

1.167 

0.060 

157. 

0.30 

0. 943 

0. 76b 

21. 7 

1 . 1 30 

26.7 

656. 

26o7 . 

0.2458 

1.333 

0. Of 6 

156. 

1.0 

0.951 

0.763 

30.1 

1-274 

37 . 4 

65i. 

3)33. 

0.1954 

1.667 

0.0’3 

1 5^. 

1.145 

0.953 

0.770 

36.8 

1.92. 

45.3 

649. 

3S16. 

0.1700 

1.903 

0.087 

154. 

1.5 

0.94 3 

0. T 64 

38.6 

2.007 

47. 3 

662- 

3S4-. 

0.1700 

2.500 

0.119 
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0.70 

0.7896 

22.44 

0.6308 

37.36 

0.7911 

0.950 
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0.75 
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above the backstreaming limit if the gap is maintained at the minimum 
value. The maximum beam current in this range is obtained by increas- 
ing the gap (at maximum electric field) until the total voltage, with 
the calculated net voltage, corresponds to the maximum permissible 5. 
This region also corresponds to constant thrust (see Figures A-3(c) 
and A-10(c)). In a space-charge-limited condition, the electric field 
force at the accelerator corresponds to the time-rat e-of-change of 
charged particle momentum between the emitter and the accelerator. With 
constant electric field, then, the thrust force at the acclerator is 
constant. With R also constant, a constant fraction of this thrust 
appear^ as thrust after deceleration. 

Although the number of grids is not indicated, a three-grid ion 
optics assembly is assumed for any R below 0.7. In the R range from 
0.7 to 0.8, either two- or thT^e-gt id optics could be used. 

The thrust-to-power ratio (see 7 igures A-3(d) and A-iO(d)) gener- 
ally rises as specific impulse is reduced, as one would expect from the 
reduced kinetic energy requirements of lower specific imp lses. What 
is less obvious is the leveling off of thrust-to-power ratio at low 
specific impulses. For a constant discharge loss and a constant pro- 
pellant utilization, the thrust-to-power ratio will actually show a 
maximum at a net voltage equal to the discharge loss in eV/ion (power 
efficiency of 0.5). The optimization of propellant utilization used 
here prevents this maximum from being evident. As the specific impulse 
is decreased below the value where this maximum might be expected, the 
optimization procedure selects increasingly lower propellant utiliza- 
tions. For example, argon at 1500 sec gives an optimum propellant util- 
ization of 0.454, while xenon at 1000 sec gives 0.613. The thrust-to- 
power ratio rises slightly at these low specific impulses because the 
low utilizations are obtained by using very shallow discharge chambers, 
which, in turn, give very low discharge losses. For argon at 1500 sec, 
the 0. 5-m-diameter chamber has a depth of only 2 cm, while for xenon 
at 100C sec, the depth is only 1.1 cm. 
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Optimum operation at very low specific impulses is, in effect, 
obtained by operating at near optimum electrical values and throwing 
away increasing fractions of un-ionized propellant as specific impulse 
is decreased. Such operation will generally not be of interest in mis- 
sion analysis: if it should appear to be preferred, it usually indi- 

cates that a stage mismatch exists. That is, overall performance would 
be increased by using a smaller mass for the electric propulsion stage 
and operating the thruster at a high specific impulse. Such a mismatch 
may exist in any case where a thruster efficiency less than about 0.5 
is indicated as being desirable. 

Sensitivity analyses are shown for variations of propellant utili- 
zation about the optimum values in Figures A-4 through A-6 for argon 
and Figures A-ll through A-13 for xenon. In the low specific impulse 
range (R = 0.17 at optimum utilization), shown in Figures A-4 and A-ll, 
relatively large departures from optimum propellant utilization cause 
small decreases from maximum thruster efficiency, not thrust, as clearly 
indicated by Figures A-4(b) and A-ll(b). Small decreases in thruster 
efficiency from maximum can be exchanged for relatively large increases 
in thrust. These thrust increases correspond to operating at higher 
net and total voltages, so that the increased thrusts result from the 
increased current capacities of the ion optics at the higher vol ages. 

At a sufficiently low propellant utilization, the total voltage can 
reach the maximum electric field value and R must increase above 0.17. 
Although shown only for xenon, the same change in R occurs at lower pro- 
pellant utilization than covered in Figure A-4 for argon. 

the effect of varying propellant utilization away from optimum has 
a more pronounced effect on thruster efficiency at higher specific 
impulses, as shown in Figures A-5, A-6, A-12, and A-13. Decreasing 
propellant utilization from the optimum value increases thrust slightly 
in the intermediate specific impulse range (Figures A-5 and A-12), and 
not all in the high specific impulse range (Figures A-6 and .-13). 

The effects of varying thruster diameter are show, in Figures A- 7 
through A- 9 and Figures A-14 through A-16. All data shown in these 
t Igures were optimized for maximum thruster efficiency at each thruster 
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diameter. As with the effect of propellant utilization variations, the 
nature of diameter effects depend on the specific impulse range under 
consideration. 

At the lowest specific impulse. Figures A-7 and A-14, the grid gap 
remains at the 0.95 mm value until a beam diameter of 0.57 m is reached. 
Above this diameter, the span-to-gap limit required that the gap be 
increased proportionally with the diameter. The effect on thrust can 
best be understood by first ignoring the effect of propellant utiliza- 
tion changes. For beam diameters less than 0.57 m, the grid gap is 
fixed at 0.95 mm, and the thrust per unit area should be constant at a 
constant specific impulse. At beam diameters larger than 0.57 m, all 
ion optics dimensions are increased in proportion to the beam diameter. 
At constant propellant utilization, the voltages would again be constant 
A classic space-charge-flow calculation has shown that beam current 
(and thrust) are independent of ion optics size when all dimensions are 
changed in proportion with voltages held constant. 

We will now include the additional effects of changing propellant 
utilization. As beam diameter is decreased from 0.57 m, the volume- 
to-area ratio, the primary electron region, VWA^, tends to decrease. 
Although (A^) also decreases, beam current decreases about the same 
amount, and hence the net effect on propellant utilization is negligible 

This tendency to decrease i /A can be partiallv offset bv an increase 

P P ' 

in discharge chamber 2 / d, . But, because an increase in t /d, also 

o b o b 

results in an increase in discharge losses, the increase in C^/d^ must 
be a compromise between maintaining propellant utilization and increas- 
ing discharge losses. This compromise results in a net decrease in 
optimized thruster efficiency. There is also a net loss in propellant 
utilization as beam diameter is decreased below 0.57 m. This loss result 
in increases in net and total voltages; these increases are the cause 
of the increase in the thrust-to-beam area ratio as beam diameter is 
increased. 

When beam diameter is increased above 0.57 m, the ratio y /A tends 

P P 

to increase. But this is more than offset by the increase in A , which 

o 
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increases with the beam diameter squared. But with beam current tending 
to remain constant as beam diameter is increased, the propellant utiliza- 
tion must suffer a net decrease. This net decrease, in turn, results 
in a decrease in optimized thruster efficiency as beam diameter is 
increased above 0.57 m. 

Maximum thruster efficiency (all beam diameters) in the low specific 
impulse range thus corresponds to the maximum beam diameter for which the 
minimum grid gap can be maintained. The thruster mass has not been dis- 
cussed, but it normally increases somewhat less rapidly than beam area. 

It appears that 'iiaximum thrust-to-weight ratio would also be found at a 
beam diameter of 0.57 m. Further, the mass penalty for going above 
0.57 m appears to be substantial. Large thrust applications in this 
specific impulse range thus strongly indicate the use of the maximum 
beam diameter consistent with the minimum grid gap. 

The tradeoffs are slightly different for beam diameter variations 
in the intermediate specific impulse range (Figures A-8 and A-15). 

Above a beam diameter of 0.57 m, the grid gap again increases above the 
0.95 mm value. However, the freedom to decrease R as the grid gap is 
increased more than offsets any thruster efficiency loss that might be 
expected from the grid gap increase. The efficiency tradeoff shifts 
when R is decreased to 0.17 since further diameter increases result in 
a net loss in optimized thruster efficiency as well as a rapid loss in 
the thrust-to-beam area ratio. 

The variation in thrust-to-mass itio is probably not great in the 

diameter range with l >0.95 mm and R > 0.17. The need for both maxi- 

8 

mum thruster efficiency and minimum parts count (for increased reliabil- 
ity and reduced cost) would thus indicate a beam diameter for large 
thrust systems such that an R of 0.17 is just reached. That is, the beam 
diameter is such that the grid gap (from the span-to-gap limit) corres- 
ponds simultaneously to maximum electric field and to an R of 0.17. 

Except or having a region with R = 0.8, the qualitative effects of 
varying beam diameter in the high specific impulse range (Figures A-9 
and A-16) are similar to the effects described above. Again, the 
optimum beam diameter appears to correspond to maximum electric field 
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and an R of 0.17 at a grid gap given by the span-to-gap limit. The 
optimum beam diameter at 10000 sec, however, is quite large — 3.90 m 
for argon and 11.76 m for xenon. The power levels are even larger for 
this optimum beam diameter — 1.26 MW for argon and 10.31 MW for xenon. 

B. CONCLUSIONS 

Much of the significance of tne large thruster analysis presented 
here lies in the specific performance values obtained. These values 
are adequately presented in the various figures and tables that have 
been included and need not be repeated here. There are other signifi- 
cant results, however, that may not be as obvious. 

One result was the extremely shallow depths obtained for optimum 
discharge chambers. At 5000 sec, fc example, the optimum depths for 
0. 5-m-diameter chambers were less than 0.1 m for both argon and xenon. 

At lower specific impulses, the optimum chamber depths were even less 
than 10 cm. The sensitivity analysis for propellant utilization indi- 
cates that small increases in depth from the optimum value will result 
in even smaller losses in efficiency. Even so, the optimum depths are 
well below what might be expected without detailed analysis. Only at 
specific impulses approaching 10,000 sec do the length-to-diameter ratios 
of the discharge chambers approach values normally expected. 

The distribution of electrons to a chamber with an & /d, of 0.2, 

P b 

or less, is another aspect that deserves mention. Obtaining a uniform 
distribution of primaries from a single baffle annulus is probably 
beyond present technology. 

The effect of beam diameter on thruster performance was examined 
at low, intermediate, and high specific impulses. The optimum choice 
of beam diameter for very large systems was indicated for these three 
ranges. The optimum choice of beam diameter can also be described in 
terms of simultaneously satisfying multiple limits. These limits are 
minimum net-to-total voltage ratio, maximum span-to-gap ratio, and 
either minimum gap or maximum electric field. Further examination of 
these limits should obviously be considered in future technology 
studies. 
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The size, power, and efficiency limits indicated by this study 
are clearly much more restrictive at low specific impulses. Alternatives 
to "conventional" discharge chambers and ion optics could offer sub- 
stantial advantages at these low specific impulses, and should therefore 
be considered. 
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APPENDIX B 


DETERMINATION OF MULTIPLY CHARGED ION EFFECTS 

Operation of the gas thruster at a specified thrust level and 
specific impulse defines the propellant flow rate to the thruster. If 
the production and loss rates of both ionic species (singly and doubly 
charged) are required to balance and the loss rate of propellant is required 
to be equal to its supply rate, then the neutral to singly ionized atom 
density ratio is gxven by 


+ 9 


(B.l) 


where the quantities n Q and n + represent, respectively, the densities of 
neutrals and singly charged ions, and the Z's represent the primary 
electron rate factors for transition from the subscripted to the super- 
scripted states. The effect of Maxwellian electrons has been neglected 
in this analysis because their effect is calculated to be small. The 
plasma parameter C defined in Eq. B.l is given by 



(B.2) 


where v+ is the Bohm velocity for singly charged ions, v is the thermal 
velocity for neutral atoms, F + is the uniformity factor for singly 
charged ions, A is the effective open area of the screen grid to ions, 
A a is the effective open area of the grid system to neutrals, and is 
the ionization efficiency defined by 


n 


i 



where F + and are the fluxes of singly and doubly charged ions that 
have left the discharge chamber through the ion beam, and r is the 
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flux of neutral atoms into the discharge chamber. For xenon propellant, 
a neutral atom temperature of about 800°K, a Maxwellian electron tem- 
perature of about 4 eV, and effective open-area fractions for the screen 
and accelerator grids of 0.7 and 0.2, respectively, Eq. B.2 becomes 


C 



(B. 3) 


The propellant utilization is left as a parameter to facilitate the 
optimization of the total design, and the uniformity factor for singly 
charged ions also appears because of the uncertainty in the proper value 
for this parameter. This uncertainty comes about because of the smell 
amount of test data available for both the multipole configuration and 
argon and xenon propellants. 

Once the ratio of densities of neutral atoms to singly charged ions 
is known, the ratio of doubly to singly charged ion fluxes can be deter- 
mined using 




(B. 4) 


This equation is a result of equating the production and loss rates of 
doubly charged ions. To evaluate it, the assumption of approximately 
equal uniformity factors for doubly and singly charged ions has been 
made (F, . = F ). Investigating the behavior predicted by Eq. B.4 for 

i T “r 

xenon, Figure B-l shows the variation of the ion flux ratio for varia- 
tion of the parameter C, with the potential difference between the anode 

and keeper (V -V ) as a parameter and using the rate factors given in 

^ ^ B— 1 

Table B-l. This parameter enters because the rate factors appearing 

in Eq. B.l are dependent on primary electron energy, which is approxi- 
mately equal to this potential difference. 
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RATIO OF DOUBLY TO SINGLY 
CHARGED ION FLUX <1"^/!^) 



0 0.4 0.8 1.2 1.6 2.0 2.4 


PLASMA PARAMETER <C) 

Figure B-l. Variation of the ratio of doubly charged ion flux to 
singly charged ion flux as a function of the plasma 
parameter C (xenoa propellant). 
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Table B-l. Xenon Race Factors 



Multiplying 'he flux ratio determined from Figure B-l by two yields 
the ratio of doubly to singly ionized currents; dividing the flux ratio 
by Jl yields the particle density ratio for the two species. This 
figure shows that the doubly charged ion flux ratio rises rapidly as C 
is reduced to very low values. This can be realized physically by 
increasing the utilization to values very near unity, a condition that 
would correspond to operation with an excessively large volume-to-surface 
area region for the prevailing densities in the chamber. Low values of 
the plasma parameter C can al30 he realized if there are steep gradients 
in plasma properties near boundaries. This condition is reflected in an 
increase in the value of the uniformity factor F + . 

Under some operating conditions for the xenon thruster, the doubly 
ch«. . ged ion density can become comparable to the singly charged ion 
density. When this occurs, the possibility of significant triply charged 
icn production becomes a concern. The dominant reaction whereby these 
ions are formed is 


xt* + e' * xt* + 2e~ . (B.5) 

Computation of the density of the triply charged specie requires as 
input the cross-section curves fo the rea*. - ion. Figure B-2 shows 
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Figure B-2. Cros.s section for production of triply ionized xenon. 
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experimental data for this reaction obtained by Kupriyanov and Latypov. 
These data suggest that the density of triply charged ions would not be 
significant because the energies required for significant production 
(^60 eV) would probably not be available. There is, however, a question 
about the suitability of these data for the application considered here, 
both because of possible errors pointed out by the authors and because 
the cross section is essentially zero in the region immediately above 
the 32.1-V ionization potential associated with the reaction. 

To determine what would appear to be a more reasonable cross sec- 

B_3 

tion, a theoretical approach due to Gryzinski was used to calculate 
them: 


0 



-10 




, (B. 6) 


where - E/U^, E is the energy of the ionizing electron (eV) , and 
is the ionization potential (eV) for the i^ 1 electron. The only elec- 
trons that need to be considered in the above summation are the four 5p 
electrons left in the outer shell of doubly ionized xenon. Evaluating 
the above as a function of energy yields the theoretical curve shown in 
Figure B-2. 

Using the theoretical curve results ot Figure B-2, the primary 
+++ 3 

electron rate factor Z (m /sec) inducing the subject reactions can be 

I T 

determined from 



(B. 7) 


Values obtained from this equation and the theoretical curve ot 
Figure B-2 will be used in subsequent analysis because they will give 
higher values of the rate factor and thereby assure a conservative 
design. 
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Equating the production and loss rates of the triply charged xenon 
gives the ratio of triply to doubly charged ion densities (n t 1 1 /n t ( ) , 



where F . . , is the uniformity factor for triply charged ions, n is the 

III P 

primary electron density, v , , is the Bohm velocity for triply charged 

TTT 

ions, and v /A is the volume to surface area ratio of the primary elec- 
P P 

tron region of the thruster. The loss rate of doubly charged ions 
resulting from the production of triply charged ions should be accounted 
for; for the case being considered here, however, the triply charged ion 
production is sufficiently small that this correction would be insig- 
nificant. From previous work, the Bohm velocity and primary electron 
density are given by 


= /3 


(B.9) 


and 


n 


P 


A v 

_E 


+ " t 


CF. 



(B. 10) 


Combining Eos. B.8 through B. 10 and using the above numerical 
values yield 



(B. 11) 


For the case of interest here, where C “ 0.65, the rate factors are 
given by Table B-l and Figure B-2, and the density ratio is computed 
using Eq. B.ll to obtain the values in Table B-2. 
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Table B-2. Triply Charged Ion Density as a Function 
of Primary Electron Energy (V D -V R ) 


V -V 
D K 

n + /n o 


n -H-+ /n -H- 

40 

1.18 

6.53 x 10~ 15 

0.027 

35 

1.09 

1.70 x 10 -15 

0.008 

32 

1.02 

0 

0 


7013 


The ratifi of triply charged ion flux (T , , . ) to doubly charged ion flux 

• ft 

(r ) is given by 



(B. 12) 


and the ratio F , , /T. is given by 

TTT + 




(B. 13) 


Because the triply charged ion flux will be small, it can be 
neglected in all of the following analysis except that associated with 
screen grid erosion. Neglecting triply charged ions and using the 
definition of specific impulse while expressing the combined thrust of 
both singly and doubly charged ions in terms of the net accelerating 
voltage (V^) and the flux of each ionic species, one obtains the 
following expression for net accelerating voltage: 


V. = 


2 2 
I m g 
sp o a o 

2e n? 


1 + F ^ /r + 

i + /2 r ++ /r + 


(B. 14) 


Substituting into this equation the xenon atomic mass [m - (131) 

-27 ° 2 

(1.67 x 10" kg], the acceleration due to gravity (g Q * 9.8 m/sec ), 

the specific impulse (I = 3530 sec), and the electron charge 

(1.6 x 10 _19 C) yields: 
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1 + V r + 

i + Ji r^/r,. 


v b n* - 818 


(B. 15) 


For the case under consideration, the results of doing this are shown 
graphically in Figure B-3, where the left side of Eq. B.15 (net accel- 
erating voltage parameter) is plotted against the parameter C. 

The beam current which would be measured from the thruster is 
given by 
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(B. 16) 


where J q is the neutral flow rate into the thruster in units of A equiv. 
This flow rate is fixed under the constraint of operation at a specified 
r ecific impulse and thrust (T) by the expression defining specific 
impulse (I ) : 


J 

o 



sp o 


(B. 17) 


The neutral propellant flow rate determined from Eq. B.17 for the 
case of xenon propellant at a 3530-sec specific impulse and a C.5-N 
thrust le/el is 10.6 A. The beam current as determined from the results 
of Figure B-l and Eq. B.16 is given in Figure B-4 as a function of the 
parameter C, The decrease in beam- current as C increases is a result 
of corresponding reductions in the doubly to singly charged ion flux 
ratio (Figure B-l) . 

The beam current and net acceleration voltage required to effect 
operation under the specified specific impulse and thrust conditions 
can also be computed; the equations for these two quantities become 
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PLASMA PARAMETER (C> 

Figure B-4. Beam current parameter variation. 
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and 


818 
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i * r^./r + + r w /r + 
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(B. 19> 


For the levels of triply charged ions encountered here, however, this 
current and voltage are not altered significantly from the values that 
neglect triply charged ions. 


It is convenient to be able to calculate the thrust of an ion 
source from measured beam currents as net accelerating voltages. The 
ratio of the actual thrust to that calculated using these quantities is 
referred to as the thrust correction factor for doubly charged ions, a. 

It is plotted against the parameter C in Figure B-5. 

The ionization efficiencies achievable using argon as the propel~ 
lant do not cause nearly as much double ionization as when the thruster 
is operated with xenon, and formation of triply charged ions is neg- 
ligible. Design curves comparable to those shown for xenon (as functions 
of the plasma parameter C) are shown in Figures B-6 through B-9 for 
operation of a thruster on argon. 
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Figure B-5. Variation of the thrust correction factor for doubly 

charged ions, a, as a function of the plasma parameter 
C (with V -V„ as a parameter). Xenon propellant. 
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Figure B-7. Variation of the net acceleration voltage 
parameter as a function of the plasma 
parameter C. Argon propellant. 
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Figure B-8. Variation of the beam current parameter as a 
function of the plasma parameter C. Argon 
propellant . 
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Figure B-9. Variation of the thrust correction factor, a, as a 
function of the plasma parameter C. Argon 
propel’ *ut. 
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